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ABSTRACT 


To make a static or dynamic structural engineering analysis of a 
structure laminated of fiber-reinforced composite material, one must 
know beforehand the macroscopic equivalent orthotropic plate properties 
of a single layer. In the case of a monofilament composite layer, 
which consists of only a single row of parallel, equally spaced rein- 
forced filaments embedded in the matrix material, flexural and twist- 
ing stiffness analyses must be carried out in addition to the in-plane 
and thickness-shear stiffness analyses owing to the inhomogeneous pro- 
perty distribution along the thickness direction of the layer. 

In this work are carried out elastic and damping analyses result- 
ing in determinations of all of the various stiffnesses and associated 
loss tangents for the complete characterization of the elastic and damp- 
ing behavior of a monofilament composite layer. 

For the determination of the various stiffnesses, either an 
elementary mechanics-of-materials formulation or a more rigorous 
mixed-boundary-value elasticity formulation is used. The solution for 
the latter formulation is obtained by means of the boundary-point least- 
square error technique. 

Kimball -Love 11 type damping is assumed for each of the con- 
stituent materials. For determining the loss tangents associated with 
the various stiffnesses, either the viscoelastic correspondence principle 
or an energy analysis based on the appropriate elastic stress distribution 
is used. 

iii 


The numerical results obtained from this analysis compare favorably 
with some existing analytical and experimental results for various 
practical combinations of constituent materials , specifically boron- 
epoxy, boron-aluminum, and E glass-epoxy composites. Finally a 
complete set of property curves for a boron-epoxy monofilament com- 
posite for a frequency range of 50 to 2,000 Hz are presented to 
facilitate elastic and damping analyses of a structure consisting of 
multiple layers of such a composite. Design data for other composites, 
namely, boron-aluminum and E glass-epoxy, which are valid for a 
similar frequency range, are summarized in tabulated form. Complete 
computer program documentation and listing are given to facilitate 
calculations of monolayer property values for other composites. 
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SECTION I 


INTRODUCTION 

1.1 Introductory Remarks 

Engineering structural design requirements for aerospace vehicle 
structures demand a maximum of strength and stiffness at minimum weight. 

In order to fulfill these requirements, engineers have been searching 
for new and better materials. This led to the development of glass- 
fiber-reinforced plastics in the case of pressure vessels, and to the 
use of carbon or boron filamentary composites in the case of stiffness- 
critical structures such as fuselage panels. In addition to being highly 
efficient structurally (strength/weight and stiffness/weight) , the 
modern composite structures present modern structural designers with a 
unique advantage over the conventional homogeneous materials in that they 
can be designed to give different properties in different directions as 
required by the particular application. This is commonly achieved by 

ic 

laminating several layers of the monofilament composites with the fila- 
ments of each layer oriented in some prescribed direction. Before 
macroscopic properties of such laminates can be predicted, it is necessary 
to determine the macroscopic behavior of a single monofilament layer. 

In contrast with those filamentary composites with many small reinforcing 
filaments randomly distributed throughout the entire cross section, there 
is only a single row of regularly-spaced filaments embedded in the midplane 
of the layer for the monofilament composites (reference 1). 




Owing to their oriented nonhomogeneous nature, one-layer filamentary 
(composite materials behave as an orthotropic material on a macroscopic 
Ibasis. Thus, to use these advantages fully, orthotropic properties of 
;the composites must be characterized from the knowledge of constituent 
material properties and their respective geometrical configurations. 

In the case of a one-layer filamentary composite with many small parallel 
:filaments more or less randomly distributed throughout the cross section, 
!the complete characterization of the in-plane macroscopic composite 
(properties requires four independent elastic coefficients. These coef- 
ficients are: two moduli or elasticity E^, ^22’ one mo ^ u * us rigidity 

* 

G ; and one independent Poisson s ratio v. _ (reference 2) . The flexural 
66 12 

and twisting stiffnesses are then related to the in-plane moduli by the 
formula (reference 3) , 

h/2 

(D H ,D 12 ,D 22 ,D 66 ) = J_ h/2 (E 11 /X ’ V 12 E 22 /X ‘ E 22 /X ’ G 66 ) z dz 

where 

X * 1 " v 12 v 21 


However, for a monofilament composite such as boron-epoxy, 
consisting of only one row of filaments whose diameter is relatively 
large in comparison with the thickness dimension, there is a large amount 


* 


In usual notations, "1" is the filament direction, "G " (i,j=l,2,3) 
denoted the shear modulus in the ij-plane, and (i,j=l,2,3) denotes 

the j-direction normal strain due to unit i-direction normal stress. 



of relatively flexible matrix material located at an appreciable distance 
from the filament axis. (See figure 1). Thus, the conventional practice 
(references 3 and A) of assuming homogeneous property distribution in 
the thickness direction will not be expected to be valid in predicting 
the macroscopic flexural and twisting stiffnesses of a monofilament com- 
posite. In view of this, flexural and torsional analyses must be car- 
ried out in addition to the in-plane analyses in order to describe the 
behavior of such a composite completely. Furthermore, owing to the 

’ k 

presence of a relatively flexible matrix material, the thickness-shear 
flexibility is expected to be significant (references 5 and 6) for the 
filamentary composites. 

In the dynamic analysis of a structure consisting of multiple 
layers of filamentary composites, the damping characteristics are just 
as important as the stiffness characteristics. The damping properties 
of a composite may be characterized in a number of ways (reference 7; 
also, see Appendix B) . 

There have been a few experimental investigations on damping character 
istics of sandwich materials and filamentary composites (references 
8-19). 

Pottinger (ref. 8) measured the temporal decay of axial vibrations 
of bars of glass fiber-epoxy and boron fiber-aluminum composites in the 
frequency range of 1 kHz to 100 kHz. 

Often referred to as transverse-shear, here the term thickness -shear 
is used, so that the term transverse can be reserved to refer to the 
direction normal to the longitudinal (filament) direction and con- 
tained in the plane of the layer. 


Schultz and Tsai (refs. 9 and 10) used the free vibration decay and 
resonant response of cantilever beams made of unidirectional and angle- 
plied glass fiber-epoxy composites in the 10 to 10,000 Hz range. 

James (ref. 11) and Bert et al. (ref. 12 and 13) used the temporal 
decay of free-free sandwich beams. The facings of the ref. 12-13 beams 
were of glass fiber-epoxy. 

Clary (ref. 14) conducted resonant response experiments of free- 
edge plates made of unidirectional boron fiber-epoxy composite material 
to study the effect of fiber orientation. 

Bert et al. (refs. 15-17) carried out resonant response measurements 
on a free-edge, circular, truncated conical shell with an aluminum 
honeycomb core and glass fiber-epoxy facings. 

Richter (ref. 18) used the rotating beam deflection technique to 
determine the damping characteristics of glass fiber-epoxy at low fre- 
quency (0.01 to 1.07 Hz.). 

Kerr and Lazan (ref. 19) made hysteresis measurements on a sand- 
wich beam with both core and facings of glass fiber-epoxy. 

Analytically, Hashin (references 20 and 21) determined complex modul 
of viscoelastic composites (particulate and filamentary) by develop- 
ing a correspondence principle which relates the effective elastic 
moduli and creep compliances of. the viscoelastic composites. However, 
his analyses are not applicable to those cases where the effective 
elastic moduli are not explicitly obtained. Since the steady-state re- 
sponse of the filamentary composite structure is of our main concern in 
this investigation, the wavelengths of the modal profiles are in general 
much greater than the filament diameter or the thickness dimensions; thus 



the Kimball-Lovell type material damping (reference 22) is assumed to 
hold. The analysis similar to reference 13 is carried out by using the 
numerically obtained stress distributions for the effective elastic 
coefficients, and the results are presented in terms of complex moduli 
for each assumed loading mode: in-plane, flexural, twisting, and so forth. 

In the following subsection 1..2, the current literature on micro- 
mechanics analyses of filamentary composites are briefly surveyed. 

In Section II, detailed elasticity and mechanics-of-materia Is analyses 
are made to obtain various macroscopic elastic properties pertinent for 
the characterization of one-layer composites. In Section III, the elastic- 
analysis results are used to calculate the damping properties of the 
composites; and in Section IV, the damping properties in terms of complex 

moduli are summarized, and some typical numerical results are compared 

/ 

with existing analytical and experimental results. 

1.2. A Brief Survey of Micromechanics Analysis of Filamentary Composite 

Materials 

Since the development of practical methods for manufacturing parallel- 
filament-reinforced layers of composite materials in the 1950's (references 
23, 24), the field of micromechanics analysis mushroomed rapidly starting 
with the pioneering analyses of Outwater (reference 23) in the United 
States and of Beer (reference 26) in Germany. By micromechanics analysis 
is meant an analysis which leads to the prediction of the macroscopic 
properties--elas tic moduli-- of the composites based only on reinforce- 
ment configurations, and the properties and volume fractions of the 



constituent materials. 


Chamis and Sendeckyj, in their recent survey of the field (ref. 4), 
listed 109 references concerned with the prediction of thermoelastic 
properties of the fibrous composites. Since a unidirectional filamentary 
composite will behave macroscopica lly as an orthotropic materials, a 
complete description of the elastic properties of the composite requires 
nine independent elastic coefficients: three Young's moduli (E^ , e 22 >E 33 ) 
three shear moduli (G^,G^^,G^^) , and three Poisson's ratios (v^ 2 » v 23 » v 3 i 
where the subscript (or subscripts) refers to the orthotropic axes of 
the composite a long which the properties are measured (see figure 1.) 

The analytical methods which have been used in previous micromechanics 
analyses may be categorized as: 

(1) Netting analysis method 

(2) Mechanics-of-materials method 

(3) Self-consistent model method 

(4) Variational method 

(5) Exact classical elasticity method 

(6) Statistical method 

(7) Discrete element (finite element) method 

(8) Semiempir ical method 

(9) Microstructural method 

In netting analysis, fibers are assumed to provide all of the 
longitudinal stiffness and the matrix material is assumed to provide the 
transverse and shear stiffnesses and the Poisson's effect. This is 


equivalent to assuming the disjointed fiber-matrix model, and thus pre- 
dicts relatively low values for E ^2 and (references 25 and 27) . 

The mechanics-of-materials methods were pioneered by Ekvall (reference 
28). In his analysis, the macro-composite properties are expressed in 
terms of the averaged stress-strain states, which, in turn, are expressed 
in terms of the constituent properties using displacement continuity 
and force equilibrium conditions at the matrix-fiber interface. In 
general, the predicted transverse and shear stiffnesses are lower than 
the experimental values. Thus, the method was later improved upon by 
using the concept of restrained matrix model in which the strain in the 
matrix parallel to the fiber is assumed to be zero (ref. 1). This method 
was later extended to account for the effect of voids in the composite 
by Greszczuk (reference 29) , and the effects of misalignment by Nosarev 
(reference 30) . 

The self-consistent model method is based on the assumption that 
the strain field of a single fiber embedded in an infinite (reference 
31) or a finite (reference 32) matrix is indistinguishable from that of 
the composite. The results of such analyses are generally accurate 
for the case of low-fiber-volume fractions only. 

The variational method is based on the energy theorems of classical 
elasticity (reference 33) in which lower and upper bounds of the layer 
properties are obtained from the theorems of complementary and potential 
energy, respectively (references 34 and 35) . The upper and lower bounds 
are very far apart for composites with high fiber-matrix-stiffness ratio 
such as boron-epoxy composites. Thus, correction factors such as contiguity 
and misalignment factor, etc. need to be brought in to obtain closer 


7 


agreement between the theoretical predictions and the experimental re- 
sults (references 34 and 36) . 

The exact classical elasticity method has many variations depending 
on the methods of solution (references 37-39) . With the exception of the 
relatively simple case of a circular fiber embedded in a circular matrix 
material, the solution cannot be obtained in a closed form and thus 
various numerical methods must be used. These are practical now with the 
aid of modern high-speed digital computers. In all cases, the problem is 
formulated with an assumption that the fibers form a regular array 
(rectangular or hexagonal) ; a solution is sought for the resulting 
mixed boundary-value problem, subjected to the usual assumption of per- 
fect bonds between the fiber-matrix interface and a set of imposed 
boundary conditions (uniform tension, shear, etc.). The elastic field 
thus obtained is then averaged over the cross section and the boundary 
to yield the desired equivalent macroscopic elastic properties of the 
composite . 

In the statistical methods, the composite is modelled by the random 
distribution of the fibers in the matrix materials. Very little success 
has been achieved by this method (ref . 20) owing to the statistical averag- 
ing process that leads to insurmountable computational difficulties. 

The discrete element method was pioneered by Foye (references 40 and 

41) for the prediction of ^22 , ^12 ,v 12’ anc * V 23" H *- s resu l ts f° r circular 
filaments in square arrays are in good agreement with those of Ekvall 
(ref. 28) and Greszczuk (ref. 29) . The method may also be applied to 
cases with nonlinear matrix behavior as well as random array arrangements 
(ref. 41); however, so far its use has been relatively limited. 


The semiempirical method which is most commonly used to date is 
due to Tsai (ref. 34). He assumed that the properties of a filamentary 
composite with non-contacting fibers may be predicted by a linear inter- 
polation between the lower and upper bounds obtained from the variational 
method. This linear interpolation was improved upon by Tsai and Halpin with 
a more refined nonlinear interpolation (ref. 3) . There are other semi- 
empirical models (ref. 4) based on the equivalent section concept, parallel 
and series connected elements, and the incorporation of certain empirical 
factors. 

The microstructural method was proposed by Bolotin (reference 42). 
Postulating that the fiber behaves as a small rod and that the distances between 
the fibers are small in comparison with the characteristic distance of the 
body and using a variational principle, he derived the displacement equilib- 
rium equation similar to those for a Cosserat medium, that is, a medium 
possessing an unsymmetric stress tensor containing couple stresses. Later, 
this microstructural model was applied by Herrmann et al.to investigate the 
transverse wave propoagation in filamentary composites (reference 43) . 

Most of the micromechanics analyses described above are based on the 

• k 

hypotheses that: (1) The fibers are regularly spaced and aligned. (2) 

The fiber and matrix materials are homogeneous and linearly elastic. (3) 

There is a complete bond at the fiber-matrix interface. (4) The com- 
posite is free of voids. (5) The composite is initially at a stress 

As a rule, these hypotheses are adhered to in the above discussed 
analyses. However, there are exceptions in which one or more hypothesis 
is relaxed, for example, perf ec t interface bond is not assumed for the 


netting analyses. 
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free state. (6) The composite behaves as a homogeneous general ortho- 
tropic material macroscopically . 

There have been a few experimental investigations on the damping 
characteristics of composites (refs. 8-19). In general, it was agreed 
upon that for small oscillations, the filamentary composites exhibit 
anisotropic, linear viscoelastic behavior. Damping properties of some 
non-metallic materials are summarized and tabulated in reference 44. 

Analytically, Hashin developed a correspondence principle which re- 
lates the effective macroscopic elastic moduli to the effective 
viscoelastic moduli (reference 45) . This correspondence principle was 
then applied to particulate and fibrous composites for the determination 
of macroscopic complex moduli of such composites (refs. 20 and 21). 
Unfortunately, this correspondence principle cannot be applied to cases 
where the effective elastic moduli are not explicitly obtained in a closed 
form. Bert etal., in their investigation of the damping in a shear- 
flexible sandwich beam (ref. 13), observed that the damping character- 
istics of composites may be related to the ratio of the dissipated energy 
per cycle to the total potential energy stored per cycle. The analysis 
is based on the assumption that the wave length of the normal modes is 
large in comparison with the thickness dimension of the beam. 
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SECTION II 


ELASTIC ANALYSES 

This section is concerned with the detailed elastic analyses lead- 
ing toward predictions of macroscopically equivalent orthotropic pro- 
perties of a monofilament composite layer. The composite layer is 
modelled by a typical repeating cross section consisting of a rec- 
tangular matrix with a centrally oriented circular-cross-section fiber. Var 
ious in-plane, flexural, twisting, and thickness-shear properties are 
obtained from the solutions of a series of mixed boundary value problems 
with approximately prescribed boundary conditions. 

2.1 Introduction and Hypotheses 

Macroscopically, a single layer of filamentary composite material 
behaves as a specially orthotropic material with respect to three mutually 
orthogonal planes: the plane normal to the fibers, the plane of the 
layer, and the plane normal to the first two planes. The intersection of 
these three planes forms three mutually orthogonal axes: the longitudinal 
(or fiber) direction, the transverse direction (normal to the fibers 
and contained in the plane of the layer), and the thickness direction 
(normal to the plane of the layer). Therefore, a complete character- 
ization of the elastic properties in three dimensions requires nine 
independent elastic coefficients (ref. 2). However, in many structural 
engineering applications of filamentary composite materials or the 
laminates of such, they are used in the form of thin panels or plates 


due to the weight considerations in such applications. In view of this, 

the thickness dimension of the composites is usually much smaller than the 

other dimensions and the radius of curvature of the structure. The three 

stress components 0^,0^, and therefore may be regarded negligibly 

small in comparison with the remaining three stress components, namely, 

* 

a 1> o 2> and . This is referred to as the generalized plane stress state. 
The constitutive equation for the generalized plane stress state is given 
by 

«li Q 12 ° 

Q 2 1 Q 2 2 0 

0 ° 2 ^66 

where the elastic coefficients are symmetric, i.e., 

Q ij - Q ji (l * ‘- 2 > 

Hence, there are only four independent elastic coefficients. Equation (1) 
may also be written in terms of the engineering moduli E^, E 2 2’ ^12* anc * 
Poisson's ratios anc ^ v 2 ^ as 


In terms of double-subscript stress notation (ref. 3. p. 16), 

a l"°ll’ a 2'°22’ 0 3 H0 33’ a 4" a 23’ CT 5"°31 ’ °6"°12 
where the x^-direction is normal to the plane of the layer. 
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°1 


' E n /k 

V 12 E 22 /X 

- 

0 


e l 

, o 2 

- 

V 21 E 11 /X 

E 22 /X 

0 

1 

e 2 

, °3, 


0 

0 

2G 12. 


he 

L 6 j 


( 2 ) 


where 

X 5 1 _ v 12 v 21 


Again, there are only four independent coefficients since the symmetry 
in the coefficient matrix, which is a consequence of the existence of the 
elastic potential (ref. 2) , demands that 

V 12 /E ll = V 21 /E 22 

It follows readily from equations (1) and (2) that, 


Q ii = E ii /X (i=1 « 2 > 
Q 12 = v 12 Q 22 



j 


Hence, for a single layer of filamentary composite with many small 
filaments more or less randomly distributed throughout the entire cross 
section, in-plane macroscopic behavior is characterized by specifying 
the four elastic coefficients Q^, Q 22 » Qj^ * and or equivalently, by 
the in-plane engineering moduli E^, ^22’ ^66’ and ° ne t * le ^ oi - sson,s 
ratios ° r V 21 * 
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Owing to the presence of the relatively flexible matrix material, 
the thickness-shear flexibility, is expected to be significant (refs. 

5 and 6) for the filamentary composites. Therefore, for complete macro- 
scopic property characterization, one will need to specify flexural, 
twisting, thickness shear stiffnesses, and flexural Poisson's ratios, 
in addition to the in-plane properties. The above-mentioned type of com- 
posite has a more or less homogeneous distribution of properties through the 
thickness. Thus, only two additional stiffnesses, namely the thickness- 
shear moduli G and G /[/[ , need to be calculated, since the flexural and 
twisting stiffnesses may be obtained by the previously stated formula, 


(D H ,D 12 ,D 22 ,D 66 ) 


r 


j 


h/2 

-h/2 


(B ll /X>V 12 E 22 /X,E 22 />l ' G 66 ) 2 dz (4) 


where ,D 10 ,D 00 , and D denote the longitudinal, Poisson, transverse, and 
twisting stiffnesses, respectively, and h = the thickness of the layer. How- 
ever, for a monofilament composite, equation (4) is not expected to hold 
because of a more predominant inhomogeneJ ty in the property distribution 
through the thickness. In view of this, for the complete characterization 
of macroscopic elastic behavior, flexural, as well as torsional, analyses 
must be carried out. 

Given constitutive properties of the constituent materials and the 
fiber-matrix geometrical configurations, determinations of these macroscopic 
equivalent orthotropic properties may be carried out in a number of ways 
as summarized in reference 4. In the subsequent analyses, approaches 
based on mechanics-of-materia 1 and classical elasticity theories are used 
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to obtain solutions which are manipulated to yield the required equivalent 
orthotropic properties. The elastic solutions obtained in this way are 
used in Section III to analyze the damping characteristics of the filament- 
ary composite materials where Kimball-Lovell type material damping (see 
Appendix B) is assumed to prevail. 

The longitudinal in-plane and flexural stiffnesses and the trans- 
verse thickness-shear stiffness are handled easily for mechanics-of- 
materials analyses. For the remainder of the elastic properties, classical 
elasticity analyses are used to formulate a series of appropriate mixed 
boundary value problems. Then these problems are solved numerically by 
means of the boundary-point-least-square method, as described in Appendix 
A, to yield the desired equivalent macroscopic properties. 

The monofilament composite material is exemplified by a repeated 
rectangular cross section consisting of a circular-cross-section fiber 
centrally located, and surrounded by matrix material as depicted in 
figure 2. 

The basic hypotheses used consistently in the subsequent analyses are 
summarized as follows: 

HI. Fibers and matrix respectively are homogeneous, linearly elastic, 
and isotropic. 

H2 . Both fibers and matrix are free of voids. 

H3. The fiber-matrix interface bonds are perfect without transitional 
region between them. 

H4. Initially, the composite is in a stress-free state and all 
thermal effects are neglected. 


15 



H5. Inertial and damping effects are neglected. 


2.2 Longitudinal In-Plane Stiffnesses 

In this subsection, two in-plane engineering moduli of elasticity, 

namely, major Young's modulus and in-plane longitudinal shear modulus 

G,s , will be discussed. The major Young's modulus E, . is estimated from 
bo i I 

the law of mechanical mixtures; whereas, the in-plane longitudinal shear 

modulus G is obtained from the result of classical theory of elasticity 
66 

analysis by Adams and Doner (reference 46). 

Major Young's Modulus E^. “ A typical repeating element of a mono- 
filament composite element is subjected to a uniform longitudinal strain 
as shown in figure 2. The longitudinal stresses induced in the fiber 
and matrix, respectively, are: 


a f E f g r 


o 

m 


- E m e i 
m 1 


(5) 


where E^ and E^ are the longitudinal Young's moduli of elasticity of the 
filament and matrix, respectively. 

The total equivalent longitudinal force P in the composite is 


P 


0 f A f 


+ a A 
m m 


( 6 ) 


where A,, and A are the cross-sectional areas of the filament and matrix, 
t m 

respectively. 

The equivalent major Young's modulus E^ of the composite is readily 


Vr 

Damping is treated separately in Section III. 
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obtained from equations (5 and 6) as 


E n ■ 

where 

A = A c + A 
f m 


The volume fractions of fiber and matrix are defined as 


( 7 ) 


( 8 ) 


V c = A /A , V = A /A 
f f mm 


(9) 


In terms of the volume fractions and V^, the major Young's modulus 
E^ is written as 


'11 


= E f V f + 


E V 
m m 


( 10 ) 


or 


E 


11 


E + (E -E ) V 
m f m f 


(ID 


since 



( 12 ) 


Since E^ > E^ in general, equation (11) shows that E^ varies linearly with 
respect to from the matrix modulus (at = 0) to the fiber modulus 
(at Vj = 1) . For a monofilament composite with a square typical element, 
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the fiber volume fraction ranges between 0 and 0.785. Therefore, the values 

of E, , /E can vary between 1 and 94.4 for E../E = 120 (borori-epoxy) . 

11 m f m r 3 

Equation (11) is usually known as the simple law of mechanical 

0 

mixtures or "law of mixtures" for brevity. This relationship is also 
valid for the cases where the filament material is transversely isotropic 
with the plane of isotropy coinciding with the cross section of the filament. 
The Young's modulus E^ in eq . (11) is then interpreted as the longitudinal 
Young's modulus of the fiber. 

Since the interaction between the constituent materials , owing to the 
difference in their Poisson's ratios, is neglected completely in this 
simplied analysis, the major in-plane Young's modulus E^ calculated from 
eq. (11) is the lower bound as demonstrated by Hill (reference 47). How- 
ever, the effects of the difference in the Poisson's ratios of the con- 
stituent materials have been shown theoretically to be minute (references 
39, 47 and 48) and confirmed experimentally. Therefore, for all practical 
purposes eq. (11) may be deemed satisfactory for the prediction of E^. 

In-plane Longitudinal Shear Modulus G ^ - Consider one quarter of 
a typical repeating element of a monofilament composite as depicted in 
figure 3. The displacement field corresponding to the applied longitudinal 
shear loading is then assumed to be of the form 

u=v=0, w=w(x,y) (13) 

with the corresponding stress components: 

t 1 = G.dwVdx , t = G^wVdy ( i=f ,m) (14) 

xz i yz i J 
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Substitution of equation (14) into the equations of equilibrium yields the 
governing partial differential equations that must be satisfied in the 
fiber and matrix regions 

[(5 2 w i /^x 2 ) + (?) 2 w*/dy 2 ) ] = 0 (i=f,m) (15) 

The boundary conditions depicted in figure 3 are 

dw m /dy = 0 along y = 0 and y = ^r 

w m = 0 along x = 0 
m - 

w = w along x = ^ir 

In solving the problem posed by equations (15) and (16) , the interfacial 
continuity conditions on the displacement and the shearing stress need 
to be considered. This leads to: 

w f = w m on C ; (17) 

and , 

G, dw*/bn= G Sw m /bn on C, (18) 

t m 1 

where n signifies the outward normal to the boundary C^. The boundary- 
value problem defined by eqs. (15-18) are then solved by the finite 
difference method (ref. 46), and the effective macroscopic shear modulus 
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of the composite is determined from 


T xz /( ^ r) (19) 

Comparisons of the shear modulus predicted by eq . (19) to those obtained 
from other analyses (refs. 35 and 38) showed that equation (19) is in 
closer agreement with the experimental values (ref. 46) . 

2.3 Transverse In-plane Stiffnesses 

Consider a typical repeating element subjected to a uniform tensile 
stress of magnitude o q in the y direction as depicted in figure 4(a). 

Because of symmetry about both coordinate axes, only one quarter of the 
repeating-element cross section needs to be considered, as shown in figure 
4(b) . 

Assuming that a state of plane strain exists in the xy (or PT|) 
plane and further that each constituent material is isotropic in this 
same plane (that is, the fibers can be transversely isotropic), one can 
formulate the problem in terms of the Airy stress function $ (reference 
49). This requires satisfaction of the following governing partial 
differential equation in the absence of body forces which vary nonlinearly 
with the spatial coordinates: 

=0, in A, (i=f,m) (20) 

4 

where 7 is the biharmonic operator, f denotes fiber, and m denotes 
matrix. 

The general solutions of equation (20) in polar coordinates (p,0) 
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are of the form (ref. 49) : 


os 

$ f (p»9) = b ^p 2 + b^p"* cos © + Y (a ^p n +b ^p n+ 2) cos n 9 (21) 

o i n n 

n=2,3, . . . 

,m, m , , , m 2 , , m 3 _ 1 m — l 

$ (p,9) = log p + b Q o + b^ o cos 9 + p cos 9 

OS 

, V / „“> n .u m n + 2 . ' m -n , . * m -n+2. 

+ / (ap+br> +ap+bp ) cos n 9 (22) 

/_i n n n r n 

n=2 ,3, . . . 

Owing to the symmetry about the x-axis, only those series terms 
which are even functions with respect to 9 are retained in equations 
(21 and 22). In the usual notations, stress, strain, and displacement 

k 

components are related to the Airy stress function i by : 

c. = p -1 (d$/dp) + p‘ 2 (3 2 $/A9 2 ) 

* 

2 2 
Og = ^ <i>/ip 

t 0 = - d/dp (p 1 d$/d9) 
c r = (1+v) E 1 C ( 1 - v) C^-v 
6 q = (1+v) E 1 [(l-v) o Q - v a r ] 

k 

, The superscripts f and m which signify fiber and matrix regions, re- 
spectively are omitted from equations (23-27) for brevity. 
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> 

> (23) 




S - (o r + 0 e )l - 0 

V r9 - 2 <1 + V) E- 1 , r9 

u^/R = (1 + v) E 1 J [(1 - ,v) O r ~v c 0 ] dp 

(25) 

u 0 /r = (1 + V) e" 1 J pC(i-v) a Q -v o ] dQ - J u r d6 

j 

where v is the Poisson's ratio. 

The rectangular components of stress and displacement components are 
then related to polar components of stresses and displacements by: 

2 2 ■» 
a = 0 cos 8 + 0. sin 0 - t . sin 2 9 
x r 0 r0 I 

I 

a - a sin 8 + 0. cos 0 + t sin 2 0 (26) 

y r 0 r0 

t =(1/2) (a -a n ) sin 2 0 + t _ cos 2 0 
xy r 0 r0 

J 

u = u cos 0 - u_ sin 0 
r 0 

v = u sin 0 + u_. cos 0 
r 0 

The unknown coefficients of the series solutions, namely, a 1 , b 1 , 

J n n 

1 i 1 i 

a n > * 1 n (i = f,m), are then determined from the symmetry conditions and 
the boundary conditions depicted in figurs 4(b), and those on C^, the 
fiber-matrix interface. 

From the symmetry of geometry and loading, it is apparent that the 
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coefficients of the odd terms in the series must vanish. 


b l ' b l 

f K f 
a = b 

n n 


a 1 =° 


m . m 
a = b 
n n 


( 28 ) 


= a 


b “‘ = 0 (n=odd) 
n 


Also, in view of the interfacial continuity conditions, the following 
relationships must be satisfied^ 


f m f m f m f m 
u r = U r * u e = u e * °r V T r0 = T r0 (29) 


Thus, coefficients b f , a *\b f ,a f ,a m may be expressed in terms of 

o n n n n 

b m and b' m . 
n n 


b f = 2 [x(l-v m )][X + (l“2v f ) ]” 1 b“ 

o m to 


a = [ - (n+1) a b + P b' ]/n 
n n n 


.f - u m 
b = ry b 
n n 


a* = [-(n+l)b“ + V b n m ]/n 


a^ m = [(l-a)b™ - (n-1) b^ m ]/n 


) (30) 


where 


X 2 G. /G 

r m 


a s X(3-4v +1) / (3-4v f +X) 

m r 


0 = 4\( 1-v ) /(X- 1) 

m 


( 31 ) 


f 
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are so chosen 


Y = [X(3-4v ) + 1]/U-1) 

m 

Finally, the remaining coefficients a m , b m , b m , and b ,m 

o o n n 

that the boundary conditions on the external boundary are satisfied at 
a discrete set of points. (See Appendix A; also, references 50-52). 

The macroscopic equivalent transverse in-plane properties may now 
be determined from the average of the displacements on the boundary. 
The plane-strain stress-strain relations of an equivalent homogeneous 
rectangular element that undergoes the same average deformation as the 
monofilament composite element are of the form (ref. 2) , 


e = 5 /E - v O /E -v O / E xy 

x xx yxyyzxzz 


o = E (v O / E + v o /E ) 
z z xz x x yz y y 


(32) 


where xy indicates that the same equation holds with the roles of x and 
y interchanged, and a superscript bar indicates the average stress, strain, 

and property values of the equivalent homogeneous orthotropic material. 

1 

Also, because of symmetry of the stiffness coefficient matrix, 


v. .E . = V, .E. 
ij J ji i 


(i,j = 


x,y,z, 


i*j) 


For the case considered, 


(33) 


o = o , a =e=0 
y o x z 


(34) 


Substitution of equations (34) into the second and the third of equations 
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(32) yields two simultaneous equations for the determination of E^ and 


e = a /E - v a / E 
y o y zy z z 


a = v O /E 
z yz O y 


From equation (35), the minor Young's modulus E ^2 and the major Poisson's 


ratio are determined as 


E = E a 2 / (E e o+a 2 ) 
y z o z y o z 


v = a /a 
zy z o 


E = E.. = E + (E, - E ) V, 
z 11 m f m f 


O - f 0 d? dT) / (u 2 6) 2 

2 z 


= dv/fty = ! (l^Sr) 


Although the major Poisson's ratio ma y be obtained from the second 
of equations (36), many experimenta 1 and analytical results (references 
28,41,48, and 53) show that the rule of mixtures may be used to pre- 
dict v, o with sufficient accuracy. In view of this, the following 


2 _> 


simplified formula for Instead of the latter of eqs. (36) will be 

used in Section 3.3 for the determination of the loss tangent associated 
with v^ 2 * 


12 


= v,V + v V 

f f mm 


(38) 


On the other hand, elementary model prediction (Reuss estimate) 
of the minor Young's modulus E ^2 results in the following expression 


E 22 * t<V E £> + ( V E „>J 


(39) 


which gives much lower values than those obtained experimentally. In fact, 
Hill (reference 54) and Paul (reference 55) showed that equation (39) 
gives the lower bound on the elastic modulus for a macroscopica 1 ly 
isotropic composite. There are some mechanics-of-ma ter ia Is analyses 
(ref. 1) which predict greater values for E 22 at low fiber volume 
fractions and smaller values for E 22 at high fiber volume fractions. 

In contrast, the classical theory of elasticity approach (ref. 56) 
predicts a value of E 22 which is consistently higher than the experi- 
mental values. 


2.4 Longitudinal and Poisson Flexural Stiffnesses 

As stated previously, for a monofilament composite layer, the 
longitudinal and Poisson flexural stiffness and cannot be cal- 
culated from the in-plane properties E ^>^22 ,v 12 anci V 21 ^ ue to t * ie 
fact that the assumption of macroscopic homogeneity of the material in 
the thickness direction is no longer valid. Therefore, in this subsection 



a relatively simple mechanics-of-materials analysis is used to obtain the ef- 
fective flexural stiffnesses. 

Longitudinal flexural stiffness D^. - An element of a monofilament com- 
posite layer subjected to a pure bending moment is shown in figure 5. Owing to 
symmetry of the loading and geometry, only one quarter of the cross 
section needs to be considered for the analysis. As a first approximation, 
assume that the Bernoulli-Euler hypothesis holds throughout the cross 
section; then the strain distribution is given as 

e = a, y (40) 

z 

where a, - longitudinal bending curvature, and z = distance from mid- 
z 

plane. 

Then, the longitudinal stress in the fiber and matrix are given by: 


o f = E f a, y , o = E or, y 
t t z m m z 


The flexural strain energy in an elemental volume one unit long 

and having cross-sectional area (a f +a ) is: 

f m 


U b = (1/2) J (0 f /E f ) da f + (1/2) J (a m /E m )da m 

’ a - a m 

f m 


or 

2U /<y 2 = (E -E )J y 2 (r 2 -y 2 )^dy + E ^ ^ry 2 dy 

o o 
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which integrates to give 


U Jot, I = (E -E )(nr 4 /32) + E (ur)V/6 (42) 

D Z I m m 

For an equivalent homogeneous orthotropic one-quarter section: 

\ 

U- /<*? = (ur) 4 6 3 /6 (43) 

O / I - 

where is the equivalent longitudinal Young's modulus for flexural 

loading. 

Equating the right-hand sides of equations (42 and 43), and solving 

for E^ /E , one obtains 
11m 

E ll )/E m =1 + (^’-D(3 tt/16u 4 6 3 ) (44) 

where 

X' s E /E (45) 

r m 

Finally, the layer flexural stiffness is given by the equation 

D U ■ E U /<, - v 12 v 21 ) f ^ 

-u6r 

or 

D u = E^ /(l-v 12 v 21 )[2(u6r) 3 /3] (46) 

Recalling that the in-plane Young's modulus E^ is estimated from the 
law of mechanical mixtures (see Section 2.2): 
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( 47 ) 


E n /E m = 1 + (X'-l) V f = 1+(\ , -1)tt/(4 U 2 6) 
the ratio E^^/E^ is expressed as 

E ll )/E ll< 1+ (X'-l) (3 tt/16u 4 6 3 > 3/C 1+(X ' -l)rr/^a 2 6] (48) 

In general, X',u, and 6 are greater than 1; thus, the ratio in 
equation (49) gives a value less than 1. This shows that for a monofil- 
ament composite, the longitudinal flexural stiffness as estimated, 
eq. (4), by using the in-plane Young's modulus will be unconservative. 

Poisson flexural stiffness . - Due to the nonuniform distribution 
through the thickness, the transverse Young's modulus E 22 and the Poisson's 
ratios and are dependent locally on the thickness coordinate of the 
composite. To carry out the integration indicated in eq. (4), ^ 22 ’ v 12 ’ and 
must be expressed explicitly in terms of the thickness coordinate. In 
the following analysis, a typical monofilament composite element is 
considered to be made from the thin dy strip depicted in figure 5(b) . 

k 

Thus, on assuming eqs. (38 and 39) to hold for the estimates of E 22 and 
of the elemental strip, one obtains the following expressions: 

For ChSy^r, (or CKTIsl) 

E 22 (T ° = E m UX */IUX -(X'-l) (1-T1 2 )*] 

& 

The Reuss estimate, eq. (39), gives excellent results for com- 
puting the effective transverse Young's modulus for a thin strip, fig. 

5(b), consisting of discrete rectangular aggregates of constituent materials. 
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J 


E n ( n ) = EjU+a'-i) d-n 2 )^]/u 

where X 1 is the moduli ratio as defined in eq. (46) and "H is the normalized 
thickness coordinate (T)=y/r) . 

In view of eqs. (49), the Poisson flexural stiffness for the monofil- 
ament composite element is calculated by the following equation: 


D 12 /D 12 = { J F(1 ° dT l + V m [(ii6) 3 -l]/[3(l-v^)]j[3(l-v^)/v m ]/(u6) 3 (50) 


where is the Poisson flexural stiffness of a homogeneous element of 

the same size and consisting entirely of matrix material, and the function 
F(T0 is defined as follows: 


D® = (2/3)E m <U6r) 3 v /(1-v 2 ) 

1 Z m m m 

F(TD - T1 2 [v 12 ODE 22 (Tl)/E ni ]/[l-v 2 2 (Ti)E 22 (lD/E 11 (n)] 


(51) 


2.5 Transverse Flexural Stiffness 


The transverse flexural stiffness is obtained in a way similar to that 

discussed for transverse tension (Section 2.3). Here, two edges of a typical 

element are subjected to linearly varying stress distributions that are 

equivalent to pure bending moments. (See figures 6(a) and (b) .) 

In view of the antisymmetry condition, the Airy stress functions in 

polar coordinates are of the form 

ao 

4 f (p,9) = b^p 3 cos 9 + £ (a*p n +b^ p n+2 ) cos n 9 (52) 

n= 3,5 .... 

in tn 3 * m L 

$ (p,9) = b^o cos 9 + a^ p cos 9 (Cont'd. on next page) 
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i 


, m , . m n+2 , 'm -n , 'm -n+2. 

(a+bp + a p +b p )cos n 9 (53) 

nn K n n 


n=l , 3 , . . . 


Equations (52 and 53) satisfy the biharmonic equations, equations (20) 
(Section 2.3). In view of equations (29) (Section 2.3), that specify 
continuity of displacements and stresses at the fiber-matrix interface, 
those coefficients defining the series solutions may be readily inter- 
related as: 


a f - [-(n+l)* b m + 0 b' m ]/n 
n n n 


b F = o' b m 
n n 

a m = [ - (n+1) b m + v b’ m ]/n 
n n n 

a^= [(l-a)b n m -(n-l)b’ m ]/n 


(54) 


where a, 3, and > are as defined previously in equations (31). 

In view of equations (54) , the matrix-region Airy stress function 
may be written as 


.m . m 3 ’ m - 1 

$ (p,9) = b^ p cos 9 + a^ o cos 9 

00 

+ y C F n (p) b™ + F n (p) b n m ] COS n 9 (55) 

n=l , 3 , . . . 


where 
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F (p) a [-(n+1) p n + np n+2 + (1-a) p n+ ^]/n 
n 

F n (p) = LYP n - (n-l)p n + np~ n+2 ]/n 

The coefficients b™, a™, b™, and b m (n=3,5,...) are then deter- 

1 l n n 

mined from the boundary conditions depicted in figure 6(b): 

o = t =0 on ? = u 

x xy 3 

a = v = 0 on 5 = 0 

x 

t = v = 0 on T1 = 0 

xy 

V= °o 5/u ’ T xy = ° on ^ = U6 

Since the second and third of equations (56) are identically satisfied by 
equation (55) , these coefficients are determined from the first and fourth 
of eqs. (56) in the sense of least square error at a discrete set of 
points on the outer boundary § = u and "n = (see Appendix A). 

The macroscopic equivalent flexural stiffness T >22 ° f a typical element 
may now be calculated from the applied bending moment M, and average weighted 
curvature ^/dy on the = U6 edge. 

D 22 = M/[95/9y] ynj6r (57) 

where 

M = 2 (ur) 2 0 /3 
o 
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(58) 


9 = (2r /I) 


f 5 tv 


] T*U6^ 


3 

I - 2\1 r /3 


1 


Finally, substitution of equation (58) into equation (57) yields the 
desired result: 


D 22 ■ < 2 { f 5 Eav/»T!3 dl}' 1 (59) 

o 

2.6 Twisting Stiffness 

To determine the macroscopic equivalent twisting stiffness D,, of a 

Ou 

one-layer monofilament composite, a rectangular cross section consisting of 
N repeating typical elements is considered. (See figure 7.) Because of the 
symmetry condition, only the quarter of the cross section, which lies in the 
first quadrant of the xy plane, needs to be considered. 

In the usual notation, u,v, and w are the displacements; and a is the 
angle of twist per unit ^-length. Then, for a small angle of twist o', we 
have (reference 57) : 

A t 

u = -oyz, v = axz, w = orcp (x,y) in R (i=f,m) 

where rp*(x,y) are the z-component displacement functions to be determined 
from the equilibrium and boundary conditions. The strain and stress 
components are readily found to be: 

33 



e =e = e = e - 0 
x y z xy 


e xz = 35 a -y]> e yz 1= ^ O' K^T 1 /By)+x], (i=f,m) 


and 



T xz = a ^[(^/frO-y] , T yz = a G^[(dcp 1 /<'y) + x] , (i=f,m) 


> ( 60 ) 


On application of the equations of equilibrium, one finds that the 
first and the second of these equations are identically satisfied and the 
third equation gives 


“2 i 
V ¥=0 


in (i = f,m) 


( 61 ) 


where 


'2 2 2 2 2 
V = (a /dx ) + (d /dy ) 


( 62 ) 


The stress equilibrium condition across the fiber-matrix interface 
requires the satifaction of 


X (dcp /dn) = 


(dcpVdn) 


+(X-l)[y (dx/dn)- x(dy/dn) ] on C. 


( 63 ) 


where n is the outward normal to the boundary C^. 
Displacement continuity requires that: 
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(64) 





on 


or, equivalently, 

dcp^/ds = dcp m /ds on Cj (65) 

where s is the arc length along the boundary C^. Similarly, on 
the vertical inter-element boundary C£ (see figure 7) , 

[acpVdSlj = [^ m /?>§]j +1 > a nd Cdcp 01 /d'Hlj = [dcp m /BTl] j+1 on c 2 

( 66 ) 

where j and j+1 refer to the jth and (j+l)-th repeating elements. 

Finally, the stress-free condition on the outer boundary (see figure 7) 
leads to 

dcp m /dn = [y (dx/dn) - x (dy/dn) J on C 3 (67) 

The aforementioned torsion problem, equations (61-67) , may also be formulated 
alternatively in terms of the complex conjugate Y 1 that satisfies the 
Riemann-Cauchy equations 

dcpVdx = di ;i /dy , dcpVdy =-SY i /Bx (68) 

It is readily shown, in view of equation (68), that the alternative 
formulation leads to: 
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( 69 ) 


a 2 • 

V f 1 = 0 , in (i=£,m) 

" f ~m 2 2 

XT = T + \ (X-l) (x + y ) + constant on (70) 

+> r * 

d'T/dn = dT m /dn on (71) 

[dT m / 9 T|, = [sT m /aTi, &Y m /ds ] j+1 on c 2 (72) 

y m = (1/2) (x +y ) + constant on (73) 

For later convenience in numerical analysis, these equations are 
normalized by introducing the following notations: 

1 = x/r, T\ = y/r, cp = cp/ r 2 , T = T/r 2 , V 2 = (d 2 /d^ 2 )+(d 2 /dTl 2 ) 

(74) 

Then, the governing differential equations and the boundary conditions 
become: 

V f =0 , V cp L =0 in (i = f,ra) (75) 

XT f =Y m + ^(X-l) (gW) 

X (dcp f /dn) = (dcp"/dn) + (X-l) [*n (d§/dn) (dTl/dn) ] 1 

/ on C x (76) 

df^/dn=dT m /dn , dcp^/ds = dcp m /ds J 

(d/d5 , = (d/d^, d/dT}) [y m ,cp m ]. +1 on C 2 (77) 
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= h ( F , 2 + ti 2 ) 

dcp"/dn = iTl (d^/dn) - § (dT]/dn) ] 
and the stress components are: 


i on C 3 (78) 


t 1 = a r [(dcpVdD-TVJ = a G i r [(d'i'Vd'T)) -T]] 


T yz 1 = a G i r [(ScpVdT^+O - -a G. r [ (SY 1 /8§) -|] 


> (i=f,m) (79) 


The solution to the problem is obtained by assuming series solutions 
for each of the fiber and matrix regions in all N elements that satisfy the 
governing partial differential equation exactly in their respective 
regions. For example, in terms of polar coordinates with the origin at 
the center of the N-th fiber, the solution for the fiber and matrix 
regions takes the form: 


<n) 


■r * l 

k=l ,2 , . . 


(n) _ 

(n) cos k 9 (n) 


a k ' P 


) (80) 


..tn 


(n) 


.(n) , V r. (n) k , ,(n) -k . Q 

b o *L [b k ’(n) tl .t '(.) lM “ 8 (.) 

k=l,2, 


It is interesting to note the fiber-region solution coefficients a^ 
may be readily expressed in terms of matrix-region solution coefficients 
b^ by the use of boundary conditions on C^, equation (76). 
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= [b + (X-l)(l-M 2 )/2]/x 

o o 

= 2b^/(\+l) + X^d 

a k = 2b k /(x+1) ( k=2 , 3 , . . . ) 

b -i = x i (b r d > 

b -k = X l b k (k=2,3,...) 

where 

\ l = (X-1)/(X+1) , d = 2 (n-1) u (82) 

In view of equations (81), only the coefficients need to be determined 
from the boundary conditions on and C^> namely, equations (77 and 78) , 
since a^ and b k may be readily calculated from b^ by the use of (81) . On 
writing b ^ in terms of b^, one has 

90 

= b^ - (X^d)p X cos 9 ^(pk+X^p k) cos k Q (83) 

k=l ,2 , . . . 

The problem is now reduced to the determination of the coefficients b^ by 
the satisfaction of equations (77 and 78) at a discrete set of points on 
C£ and in the sense of least square error (see Appendix A) . 

After determining b^> one can calculate readily the stress components 
from equation (79). Finally the torsional stiffness D,, is obtained from 

. DO 

r- / \ 

Superscripts ' are omitted here for brevity. 
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the following equation in a closed form. (See Appendix C for details.) 


k (5 2 +-n 2 ) ] dF d T\ ] 
(84) 

2.7 Longitudinal Thi ckness -Shear Stiffness 

Owing to the presence of relatively flexible matrix material, the 
effect of thickness-shear deformation needs to be considered in the struc- 
tural application of monofilament composite materials. As the thickness- 
shear stresses are not distributed uniformly across the cross section, the 
effective thickness-shear strain needs to be known for the calculation 
of effective thickness-shear stiffnesses. It is a commonly accepted practice 
to relate the effective shear strain to the average shear strain by means 
of a correction factor K which is usually referred to as the shear coef- 
ficient (reference 58) . 


D,, = (-y t + x t ) dx dy/a 

66 J J xz yz 


A c +A 
f m 


■ 2 c „ If CV * <? 2 + -" 2 )l dE - ^ + IJ ft. - 


m 


K'= Y /y 

average effective 


(85) 


Since the thickness-shear strain distribution is dependent on the shape of 
the cross section on which the thickness-shear stress acts, it is also re- 
ferred to as shear shape factor. 

In 1921, Timoshenko derived a theory of flexural beam vibration in 
which the effects of rotatory inertia as well as that of thickness shear 


were taken into account (reference 59). The shear coefficient K was 
defined as the ratio between the average shear strain and that at the 
midplane. This yielded a value of 2/3 for a homogeneous, rectangular- 
section beam. However, this value of the shear coefficient gave poor 
agreement with experimental results. In view of this, numerous 
attempts were made to obtain a better value for K which would be in 
closer agreement with experimental results. Based on the high-fre- 
quency mode of beam vibration, Mindlin and Deresiewicz (reference 60) 
obtained a value of 0.822 for a rectangular cross section; whereas, 
based on the static mode, Roark (reference 61) gave a value of K of 
5/6. Recently, Cowper (reference 62) used a different static approach 
to derive a formula for K which is in good agreement with those obtained 
by other investigators. This latter approach is deemed to be satisfactory 
for long-wavelength, low-frequency deformations such as those encountered 
in the vibration of monofilament composites. 

Therefore, in this subsection, Cowper's analysis (ref. 62) is 
extended to the nonhomogeneous case consisting of a typical repeating 
element of a monofilament composite to obtain the effective thickness- 
shear strain. The thickness-shear stiffness, which is defined as the re- 
sultant shear force divided by the effective shear strain, can then be ob- 
tained as a direct consequence of the analysis. 

First, the thickness-shear distribution is obtained from the analysi 
of a tip-loaded monofilament cantilever beam shown in figure 8(a). 

Denoting the displacement components as u,v, and w, one defines the 
mean displacements of the cross section and mean angle of rotation of the 
cross section q) by the following equations: 


u 



u dx dy 




W - (1/A) JJ w dx dy 

A 

cp = (1/I y ) JJ xw dx dy 

A 


> ( 86 ) 




where A denotes the entire cross-sectional area, and I is the moment 

y 

of inertia of the cross-sectional area with respect to the y axis. Then 
the actual displacements of a point on the cross section are written as 


u = U + u , w = W + xCp+w (87) 

where u and w are the residual displacements which are equal to the de- 
viations of the actual displacements from the weighted mean displacements. 
In view of the definitions, equations (86), 


II 


u dx dy 


-fl 


w dx dy 


-JJ 


xw dx dy = 0 


( 88 ) 


The stress-strain relation 


T 

xz 


/G 


+ w. 


(89) 


where a comma represents partial differentiation with respect to the 
spatial variable that follows it, may now be written as 


W 


, 2 + cp = <t x2 /G) - 


W, 


(90) 



where the shear modulus G is to be interpreted as that of the fiber 
or matrix material depending on whether the material point is located 
in the fiber or the matrix region of the cross section. Finally, the 
integration of equation (90) yields the desired kinematic relations among 
the mean values of the midplane slope W,^; flexural slope cp; and the 
effective shear strain v 

ef f 

W, z + cp = (1/A) JJ [(t xz /G) - w, x ] dx dy = y gff (91) 

A 

For a tip loading, the shear force Q is uniform along the entire length 
of the beam, hence, the first part of the integral in equation (91) is 
evaluated as 


If <\z /a) ■ (< W + <V G „> (92) 

A 

where Q„ and O represent the respective shear forces that act on the 
f m 

fiber and matrix regions and are related to the total shear force Q by 

Q-Q £ + Qit (93) 

In view of eq. (87), the remaining term in the integral of eq. (91) is 

J*J w >x dx dy = jj( w > x " ^P) dx dy (94) 

A A 

where cp is defined in the third of eqs. (86). Combining eqs. (91-94), 
one obtains the effective thickness-shear strain expression, 
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( 95 ) 


T ef f 


= W, + cp = (1/A) 


L(Qf /G f)+(Q m /G m ) 


-JJ. 


(w, x - cp)dx dy] 


where cp and hence w, the displacement field, needs to be known in terms 
of the parameters defining the geometrical configuration and the con- 
stituent properties before the integration can be carried out. Assuming 
that the deformation of the cross section can be approximated by that of 
a tip-loaded cantilever, and allowing the constituent materials to be 
transversely isotropic with the plane of isotropy normal to the fiber 
axis, one can readily formulate the problem using the Saint-Venant semi- 
inverse method (refs. 57 and 63). 

First, displacement components are assumed to be: 



= B v i (-t-z) (x 2 -y 2 ) + /, z 2 - (1/6) z 3 ] 


= B v i (A-z) xy (i=f,m) > (96) 

= -B { x (Iz - h z 2 ) + x 1 + C(E i /G.)-2v i ](^xy 2 )} 


where B is a constant to be determined from the boundary conditions; 

is the Young's modulus in the fiber direction; is the shear 
modulus in the vertical plane parallel to the fiber axis; is to be 
interpreted as the Poisson's ratios v 3 ^ =v 32 » anc ^ X 1 = X* ( x >y) ( i =f ,m) 
are functions to be so chosen as to satisfy equilibrium conditions 
in their respective regions. 

The stress components are readily calculated from equation (96) 
to be : 
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= - BG (x* + (1/2) v. x 2 [(E /G.)-3v.](%y 2 )] 

IX JL ll_l 

( i=f ,m) V (97) 

T yz = ‘ BG J *’y + f (E i /G i >- v i 3 x y] 

°z = " BE i x 

If the constituent materials are isotropic, E^, G^, and v. respectively 
of the constituent materials are related by 



which must be satisfied in the respective regions A r and A . 

r m 

Consideration of displacement continuity at the fiber-matrix inter- 
face requires that: 

f m f mf m , 1 

u — u , v = v , w = w on Cj (100) 

Apparently, the first and the second of equations (100) cannot be satisfied 
unless the two Poisson's ratios are equal. However, the interaction 
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between the constituent materials due to the differences in the Poisson's 


ratios has only weak effects as evidenced by many theoretical analyses 
(refs. 39, 47, and 48). Hence, it will be assumed for the subsequent 
analysis that 


v, = V = V (101) 

i m 

With this assumption, equation (101), the first two of equations (100) are 
identically satisfied and the third leads to 

X f =x m on Cj (102) 

Continuity of surface traction at the interface requires that 

t^ z (dx/dn) + t^ z (dy/dn) = t™ z (dx/dn) + T^(dy/dn) (103) 

where n denotes the outward normal coordinate to C^, hence, (dx/dn) 
and (dy/dn) are the direction cosines of the unit normal vector to the 
interface C^. 

In terms of the stress components defined in equations (97) , the 
condition of equilibrium at the interface, equation (103), is cast readily 
in the following form: 

G (d\ f /dn) - G (d x"7dn) 
f ^ m 

=- (G^-G^) £[% v x ^ + ( 1 ~h v) y^](dx/dn)+(2+v) xy(dy/dn)j 

on (104) 

The condition that the lateral surface is free from surface traction 
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leads to 


(d \ m /dn) = - v x^ + (1-^ v) y^](dx/dn) + (2+v) xy(dy/dn)| 


on C. 


(105) 


The solutions to the problems posed by the governing differential 
equations (99) and the boundary conditions, eqs. (102,104,105) are now ob- 
tained by assuming a pair of series solutions which are harmonic in the 
respective fiber and matrix regions. The coefficients of these series 
solutions are determined so as to satisfy the boundary conditions in the 
least-square-error sense (see Appendix A, also Section 2.6.) 

Since the solution procedure is very similar to that of Section 
2.6, it will be high-lighted by listing these equations which are pertinent 
to the solution. 

With the introduction of the following transformation, 


= x/r, Tl = y/r, \ = 


- i , 3 

“ X /r 


( i=f ,m) 


\ 


o 9 9 9 9 

V = (ft /ft r) + (* W) 

eqs. (99), 102, 104, and 105) are rewritten in the 
is convenient for numerical analysis: 


r doe) 


following form, which 


i2 i 


X = X 


= 0 

in A . ( i = f ,m) 
l 

(107) 

m 

on Cj 

(108) 


X(d X f /dn) -(d X m /dn) = -(X-l) [lhv% 2 + (l-%v)1) 2 ](d§/dn) 


on Ci 


+ (2 + v) f Tl. (dTj/dn)} 


(109) 


46 


(dx m /dn) = - [lh V £ 2 +( l~%v) T 2 ] (d|/dn) 

+ (2 + 0) 5 -n (dTl/dn)j on C 2 (110) 

The series solutions, which are harmonic in the respective regions, 
are assumed to be: 

CO 

f , V k > o 
X = a o + A a kP cos k 0 

k=l 

. (HD 

\ m = b o + + ^-k p cos k ® 

k=l 


In view of eqs . (108 and 109), which Warrant displacement con- 
tinuity and stress equilibrium conditions, the coefficients a^, a^, 
and b (k = 1,2,...) may be expressed in terms of b as: 

* K K 


a = b 
o o 

a^ = 2 b^/ (X+l) - X 1 (3+2v)/4 
a 3 = 2 b 3 /(X+l)+ ( Xj/4) 

= 2 b k /(X+l) (k = 3,4,...) 

b_ 1 = - Xj [b^ + (3+2v)/4] 

b -3 = " X 1 ^ b 3 

b_ k = - Xj t> k (k = 3,4, . . .) 

where 

x x = (X-1)/(X+1) 




\ ( 112 ) 
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Because of the anti-symmetry of the displacement component w with 
respect to the ^ axis, it may be shown readily that those coefficients 
with even subscripts are zero. The condition that the lateral surface 
is free from surface traction, eq . (110), takes the form: 


00 

I 


k b^ [pk 1 cos (k-1) 9 + Xjp k 1 cos (k+1) 9] 


> 


0 / 

- ~ (\^M) [(3+2v)p cos 2 9 - 3p cos 4 9] 

- u-l) [ % v % 2 + (l-v v)T\ 2 ] 

(0s9^tan ^6, §=U , 0ST^u6) 

00 

k b^ [-p^ 1 sin (k-1) 9 + X^p ^ 1 sin (k+1) 9] 
k-1, 3, 5 

= - (X^/4) [(3+2v) p 2 sin 2 9 - 3p ^ sin 4 9] 


\ 

f 


| 

I 


(113) 


- (X-l) (2 + v) % t\ 

(tan" 1 5<;9<n/2, 0S§<U , T)=u 5) 


Next the coefficients b^ (k « 1,3,...) are obtained according to 

the boundary-point least-square method as described in Appendix A. 

With the coefficients b. thus obtained, the constant B which 

k 

appears in eqs. (96 and 97) is calculated from the condition that the 
resultant shear force is equal to the externally applied tip load Q: 


IJ \z dx dy ■ Q 


(114) 
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Since the stress components in eqs. (97) satisfy equations of equili- 
brium in the absence of body forces, the following relation must hold: 


t = BE. x = 

ZZ,Z 1 


- t - t (i=f,m) (115) 

zx,x zy,y 


or 


T 1 + t 1 + BE. x = 0 (i = f.m) 
zx,x zy,y l 


(116) 


In view of eq. (116), eq. (114) may be written as 


q ■ 1 If. \z d * 


dy 


i=f,m i 


= Y [I [(* T 1 ) , + (y T 1 ), ] dx dy 

L J J zx' ’x zy y 

i=f,m A. 


+ B I E i JJ x2 dx d y 

i=f,m A. 

9 i 


= B 


I ‘til 

i=f,m A. 

l 


x dx dy 


s B I, 


(117) 


where A^ and A are the cross-sectional areas of the respective fiber 
f m 

and matrix regions and I is the weighted moment of inertia of the 
cross section. From eq. (117), the constant B is readily obtained as 


B = Q/l, 


( 118 ) 


where 
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x e = I E i II x2 dx dy < 119 > 

i=f,m A. 

1 

Finally, the longitudinal thickness-shear flexibility S,.^. is obtained 
from eq. (95) as 

S 55 - [Q/qi6r)]/Y eff 

= A I E ^ I x -I y^ ■ ( A / X y) IJ x (^. 1 +xy 2 )dx dy]' 1 (120) 

A 

For an equivalent homogeneous beam with a rectangular cross 
section, the longitudinal thickness-shear stiffness is given by 

S 55 = KA G 55 /(2HAr) (121) 

where the shear coefficient K is obtained as (ref. 60) 

K = 10 (1+0) /( 12+1 lO) (122) 

Thus, on equating eqs. (120) and (121), one obtains the equivalent 
shear modulus G,-^ as 

G„ = I_ (12+llv) [ 10(1+0) l" 1 [%v(I -I ) 

jj c x y 

- (A/I ) x (\ l +xy 2 ) dx dy] 1 (123) 

y JJ A 

2.8 Transverse Thickness-Shear Stiffness 

A typical repeating one-quarter cross section shown in figure 9 
is considered. As a first approximation assuming that the Bernoulli- 
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Euler hypothesis holds, we make an elementary Jourawski-type mechanics- 
of-materials analysis to obtain the shear stress and strain distri- 
butions in the cross section due to the application of shear force Q 
along the vertical boundary of the cross section. (See figure 9(a).) 
Then the strain energy is calculated and Castigl iano ' s principle is 
applied to obtain the macroscopic shearing strain v . The transverse 

thickness-shear stiffness S (/ is then obtained in terms of definite 

44 

integrals arising from the strain energy calculations. 

Assuming that the cross section remains plane, one obtains the 
flexural strain distribution as 


c x = - «. x y 


(124) 


where e denotes the flexural strain of an element located at a distance 
x 

y from the midplane, and rv >x denotes the bending curvature of the mid- 
plane. The corresponding stress distribution is 


0 L = - ,-v, y E (i=f,m) 
x x i 


(125) 


Summing the moment due to the stress distribution, one finds that 
the flexural moment acting on the cross section is 


M = 
x 


•3 -1 ? 1/2 

(2/3) E m (u6r) J [l+a’ -l)(u5) (1-0 ' ]«, x 

0^1 


(126) 


(2/3) E m (ji6r) a, x 


1 < 5*1 


where 
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( 127 ) 


= E /E , E = x/r 
t m 


From equation (126) , the flexural stiffness of a strip Ax in 
width is found to be 


°22 * \ 


f (2/3) E m (u6r) 3 [ 1+(X ' -1) (u 6) ' 3 (l-? 2 ) 3/2 ] 


(2/3) E (u6r) 
m 






(128) 


In view of equations (125-128) , the flexural stress distribution 
is now given by the expression: 


\ ■ - <M x /D 22 ) » E i 


(129) 


The thickness-shear distribution T is then obtained readily 

xy J 

from the equilibrium of the forces acting on the strip. (See figure 
9(b).) 


U6r 


U6r 


-J [° x ‘] x + J iVlxW' 5 ' ' T x, 


Ax 


(130) 


Solving for t in equation (130) and using equation (126), one obtains 
xy 


t = (A M /Ax) D I f 1 E. y dy 
xy x 22 J i 

y 


(131) 


where 


AM = [M ] -[Ml 
x x x+Ax x x 


(132) 


Since the shear force Q = A M /Ax and the integral on the right- 

- x x 

hand side of equation (131) can be evaluated, the thickness-shear 
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distribution is given by: 


For CKx<r 

t = Gv 
xy f xy 

f (Q x /2D x )E m r 2 [(X , -l)(l-5 2 )-X , 'n 2 +(u6) 2 ] (* y* (r 2 -y 2 ) % 

= l (133) 

(Q J2D ) E r 2 [ (u 6) 2 - 71 2 )]; (r 2 -y 2 ) ^<y^U6e 

x x m 

and, for r< x^ur 

t = G y =(Q /2 D ) E r 2 [(u6) 2 -T) 2 ]; 0SvKu6r (134) 

In view of equations (133 and 134), that part of the strain 
energy per unit z-length due to shearing stresses is readily cal- 
culated . 

U = (l/2)[ I T Y dx dy 

s JJ xy xy 

f m 

= [3Q 2 /(5G n p6)](u-l)+(9Q 2 /8G f )F X + (9Q x /8G m ) F 2 ( 135 ) 

where 

F 1 [ (X ' _1) (1_t \ 2) + ( U6) 2 ] 2 ( i-l) 2 ) ^ 

97 9 3/2 

- (2/3)[(X'-l)(l-T|V(u6Hx'(l-Tr) 

+ (X ,2 /5) (1-T| 2 ) 5/2 } • [(u6) 3 +(X'-l)(l-T1 2 ) 3/2 r 2 dll (136) 
F 2 = J {(8/15)(u6) 5 -(1-T1 2 )^[(u6) 4 -(2/3)(u6) 2 (1-T1 2 )+(1/5)(1-T1 2 ) 2 


Using Castigliano's theorem, one can compute the mean macroscopic 
shearing strain Y x ^ fro® equation (135) as 


by 


v xy - (^ s /SQ x )(ar) 


-1 


= (3Q x /2Ur) {[2/(5 GjjS) ] < u -l)+(3/4 ) G^Fj+O^G^fJ (138) 


The average transverse thickness-shear stiffness is given 


S., = Q /v = KAG.. 

44 x xy 44 

= (2ur G m /3)[(2/5)(U-l)/<u6) +(3/4) (F^+F^ ] _1 (139) 


where 


X = G /G 
f m 


(140) 


For a homogeneous rectangular corss section, Cowper's shear 
factor K is given by eq . (122). Thus, the equivalent transverse 
thickness-shear modulus G^, is given by 


G 44 = S 44 /(2U6r K) 


(141) 



SECTION III 


DAMPING ANALYSES 


This section is concerned with the analyses leading toward the char- 
acterization of dynamic properties of a single layer of monofilament 
composites. Many solids, whether they are metallic or otherwise, exhibit 
a damping effect ; this phenomenon is essentially due to the dissipation 
of energy associated with the deformation. It is convenient, especially 
in the case of steady-state vibrational analysis, to represent the dy- 
namic properties of the solid by the use of complex dynamic moduli . 

For example, the major Young's modulus E^ may be considered to be of 
the form 


'11 


„ R , . I 

E 11 1 E u 


(142) 


R I 

where i = /^T , E^ is the storage modulus and E^ is the loss modulus , 
since the latter is associated with that component of the strain 90° 
out of phase from the stress component which gives rise to the energy 
dissipation (reference 64). Alternatively, equation (142) may also be 
written as 
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where g is referred to as the loss tangent which is related to the 
E 11 

storage and loss moduli by the equation, 
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g E u = (E n /E n> = tan Y 

Here the term, y » *- s th e P^ ase angle by which the sinusoidally varying 

strain component lags behind the accompanying stress component on the 
phase plane, and is referred to as the loss angle. (See Figure 10.) 

In the following subsections, detailed analyses based on the results 
of Section II are made to obtain the loss tangents associated with all 
the stiffnesses necessary for the characterizing of the damping properties 
of a single layer of monofilament composite. 

3. 1 Introduction 

As previously discussed in Section 1.1, one of the main advantages 
that composite materials have over conventional homogeneous materials 
is their macroscopic anisotropic nature which gives the modern structural 
designer the ability to design stiffness and damping properties and strengths 
to fit the requirements of the application by means of lamination. In the 
dynamic analysis of a structure consisting of multiple layers of composites, 
as in the case of elastic characteristics, damping characteristics of a 
single. layer must be known beforehand in order to predict the dynamic 
behavior of the laminated structure. 

There have been numerous experimental investigations on damping 
characteristics of sandwich and filamentary composite materials (ref. 8-19). 
Analytically, Hashin (ref. 20-21) derived a correspondence principle for 
the determination of complex moduli of viscoelastic composites (particulate 
and filamentary). This correspondence principle is particularly suitable 
when the elastic moduli are of a form explicit in terms of the moduli of 
the constituent materials. However, Hashin's approach cannot be applied 
to cases where the moduli are. given in terms of numerical elastic results. 


56 



Kimball and Lovell (ref. 22) observed that many engineering materials 
undergoing sinusoidal motion exhibited energy losses which were proportional 
to the square of the amplitude and independent of the frequency. Bert 
et aL (ref. 13), in their investigation of damping in sandwich beams, 
assumed Kimball -Love 11 type damping for the facing and core materials 
and obtained good agreements between the theoretically-calculated and 
experimentally-obtained values of the logarithmic decrement for free vibra- 
tion. 

It is noted that the ratio of the total damping energy dissipated per 
cycle U d to the total potential energy stored per cycle U can be related 
to the logarithmic decrement 6 as follows (reference 65): 

6 = (1/2) ln[l-(U d /U)] (145) 

For most structural materials, where the logarithmic decrement is small, 
it may be approximated by 

6 = (l/2)(U d /U) (146) 

In turn, the loss tangent g may be related to 6 by (See Appendix B for 
details . ) 


g = (l/2n)(U d /U) (147) 

The total potential and dissipative energy per cycle of a solid under- 
going sinusoidal motion are given by * 

» - iJJkjV* <148 > 

“d - (iw 

*Repeated subscripts denote summation where i,j = 1,2,3. 
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where o_ and e denote respective stress and strain amplitudes, 

C, the damping coefficient, and V the total volume (fiber and matrix) 
d 

occupied by the composite. 

In performing the integration indicated in equations (148 and 
149) , the stress distributions are approximated by those obtained from 
the elastic analyses of Section II. Inherently, the stress distributions 
are dependent on the excitation frequency; however, for the frequency 
range of interest here (well below any micro-scale resonance), the 
stress distributions obtained from the elastic analyses of Section II 
are expected to be valid. 

For the monofilament composite element under consideration 

here, with the exceptions of in-plane longitudinal shear modulus G,, 

bo 

and the Poisson's ratio . the following general procedure is adopted 
for the calculation of loss tangent of the composite element: 

1. First, with the composite element subjected to an appropriate 

* 

loading , the strain energy is calculated for the respective fiber and 
matrix regions (U f and u" 1 ) . 

2. In view of equation (147), the damping energies for the re- 
spective fiber and matrix regions are obtained as follows: 

= 2rr g. U 1 (i=f,m) (150) 


* 

By an appropriate loading, we mean the type of loading which 
is appropriate for the property under consideration, i.e., longitudinal 
tension for E^, etc. 
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3. Then, the composite Loss tangent associated with the pro- 


perty under consideration is computed as 

g = (1/2 tt) ( i— f ,m) (151) 

i i 

The determinations of the real parts of the complex moduli; 
namely, the storage moduli; are carried out in an analogous manner to 
that of Section II by simply replacing the elastic property parameters 
with their respective counterparts in the complex moduli. 

As for the determinations of the loss tangents for the in-plane 
longitudinal shear modulus and the Poisson's ratio Hashin's 

correspondence principle is used. The equivalence of the correspondence 
principle and the general procedure described above is substantiated 
analytically for the case of the longitudinal Young's modulus and 
numerically for the case of the transverse Young's modulus ^ 22 ' 

In the subsections that follow, detailed analyses are carried 
out for the loss tangents corresponding to all of the stiffnesses 
characterizing the composite properties. 


3.2 Longitudinal In-plane Loss Tangents 

Two loss tangents, namely and ®G^ 5 > be obtained in 

this subsection. As previously discussed in Section 2.2, the major 
Young's modulus is obtained explicitly from the law of mixtures, 

whereas the longitudinal in-plane shear modulus is based on the 

results of the elasticity analysis of Adams and Doner (ref. 46) . 



Loss tangent gg^ associated with the Major Young's modulus E^ 
In terms of the notations used in Section 2.2, the total strain energy 
per cycle in the respective fiber and matrix regions are calculated as 
follows: 


U 1 = (1/2) 
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In view of equations (152 and 150) , the total strain energy per cycle 
in the composite and the damping energy per cycle are readily cal- 
culated as follows: 


u = u f + u m = ^ (AO + An) e i 
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where g, and g are the loss tangents associated with the longitudinal 
1 m 

Young's moduli of the fiber and matrix materials. The composite loss 
tangent gg^ associated with the major Young's modulus is now 

readily calculated on substituting equations (153 and 154) into equation 
(151) . 


g ~ = (g,A a + g A a ) / (a A +a A ) 

°E^ fff mmm ffmm 


(155) 


In view of equations (5 and 9), equation (155) may be written 


in terms of volume fractions V c and V and X' as follows: 

f m 



gp = (x* V g + vg)/(\’ v + V ) 

h,^ It m m t m 
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where 


X' = E R /E R 
f m 


(157) 


The composite storage modulus is obtained from eq . (10) 


as 


E R = V E R + V E R 
11 f f mm 


(158) 


It is interesting to note that identical results for the 
storage modulus and the loss tangent can also be obtained on utilizing 
the elastic-viscoelastic correspondence principle of ref. 21. Re- 
placing the elastic moduli in the right-hand side of eq . (10) with the 
corresponding complex moduli, and separating the real and the imagin- 
ary parts, one obtains 
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= (V,E R +V E R )[1 
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+ i (g,V,E R +g V E R )/(V 1: E R +V E R )] 
& f f f °m m m' f f m nr J 


(159) 


from which equations (156 and 158) readily follow. This substantiates 
the equivalence of the correspondence principle and the formula (151). 


Loss tangent associated with the in-plane longitudinal 

R 

shear modulus " As mentioned previously, the storage modulus 

is approximated by the elastic analysis of Adams and Doner (ref. 46). 

To calculate the loss tangent gG^ according to eq. (151) would require 
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lengthy numerical analyses for the determinations of the stress and 
strain distributions. To avoid this, it is surmised here that the 
results of ref. 47 may be empirically approximated by a much simpler 
formula due to Hashin and Rosen (ref. 35), which is given here in 
modified form as follows: 


G 66 /G m = °V CX(l+V f ) + (l-V f )]/[x(l- v f )+(l+V f )] 

(160) 


where Cy^ is an empirical factor brought in here so that the shear 
modulus as calculated by equation (160) coincides with that given 
in ref. 46, A is the fiber and matrix shear-moduli ratio defined in 
eq. (31), and is the fiber volume fraction. 

In view of the explicit expression of eq . (160), the loss 
tangent gG^g, * s readily obtained on replacing A by its complex counter 
part . 


A = X R (1 + ig ) 

where 


X R - (gJ/G R ) (l+g g G )/o + g 2 ) 

f m m 

g, = (gp - g r ) / (i+g r s r ) 

A 6 G G G 


(161) 
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Substitution of equations (161 and 162) into eq. (160) yields: 


G 66 /C m R ' C V f { (1-V £ 2 )Cx r2 <1-sJ)+ 1] + 2l R (l+V f 2 ) 

- 4X R 8 x g G V £ } ■ 
m 

• {[X R (1-V f ) + (1+V f )] 2 + Cx R (l-V f ) §x 3 2 } _1 (163) 

2 

g G [(l-V f 2 )[X R (l-g^ 2 )+l] + 2X R (1+V f 2 )] 
m 

+ 4X R g x V f } • 



* { (l-V f 2 )[X R (l-g 2 )+l]+2\ R (l+V f 2 )-4\ R g^g G Vf }- ] 

m 

(164) 


3.3 Transverse In-plane Loss Tangents 


Based on the result of elastic analysis and the observation made 
in section 2.3, two loss tangents &E-22 anc * Svi2 assoc l ate< ^ with the in- 
plane transverse Young's modulus E ^ anc * th e major Poisson's ratio 
will be obtained as follows. 


Loss tangent &E 22 assoc l ate ^ with the transverse in-plane Young's 

modulus . - Owing to symmetry in the loading depicted in fig. 4, shear 

stresses t and t vanish. Hence, the total strain energies per cycle 
xz yz 

for the fiber and the na trix regions, respectively, are: 


U 1 -[1/(4G )]!; [<1-V.)(CT x V)-2v.O O +2r 2 ]dx dy 

L • * 1 x y 1 x y 


(i = f ,m) 


(165) 
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where the 
5 , and 6) . 
cycle are: 


stress components 0 ^ 
In view of equations 


,a ,t are given by the equations (C-4, 
y xy a j -i > 

(150 and 165), the damping energies per 


U d = ( n S i / 2G i )j|[ (1 “ v i) (a x +0 y) " 2v i a x a y +2T xy 2 -^ dx dy (i=f > m ^ (166 > 

A. 

l 

The loss tangent &E 22 the* 1 readily calculated from eq. (151). 

The storage modulus E 22 is calculated from the first of equations (36) 
by means of appropriate storage moduli of the constituent materials 
corresponding to the particular frequency of excitation. 

Loss tangent g v ^ associated with the in-plane major Poisson's 
ratio . - Since the in-plane major Poisson's ratio is related to the 

constituent Poisson's ratios V m and by a simple law-of-mechanica 1- 
mixture formula, eq. (38), the complex Poisson's ratio i s readily 

obtained in a manner similar to that used for E^, eq . (159), by the use 
of Hashin's correspondence principle as follows: 


Vl * = ( Vf +V m V m ([ 1+1(8 
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V,v*+g V vV(V.v!+V v R )] 
V, f f V ram ff mm J 
-f m 


(167) 


where 


R „ R . „ R 
v,, = V.v f + V v 
12 f f mm 


(168) 


g = (g V.v R + g V v R )/(V v, R +V v R ) 
V, o V, ff v mm ff mm 

12 f m 


3.4 Longitudinal Flexural Loss Tangents 


In this subsection are treated damping analyses for determining 
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the loss tangents gp^ an< * g Di2 associate ^ with the longitudinal and 
Poisson flexural stiffnesses. The former is obtained according to eq. 
(151), whereas the latter is obtained by the application of Hashin's 
correspondence principle. 

Loss tangent go^ associated with the longitudinal flexural 
stiffness . - In view of eq . (42) derived previously in Section 2.4, 
the strain energies per cycle in the fiber and matrix regions are 
readily obtained as follows: 



g D n = {C(X' gf - g m )TT/ 3 2 J+ (g m 6 3 U 4/6) } 

• {rU , -l)(n/32)+ (6 3 u 4 /6)}- (170) 

The storage stiffness is readily calculated from eq . (46) with 

the help of eq. (44) by replacing the static moduli and the Poisson's 
ratios with their respective real parts of the corresponding complex 
moduli and ratios. 

Loss tangent associated with the Poisson flexural stiffness. - 

According to Hashin's correspondence principle, the complex Poisson 
flexural stiffness may be calculated from equation (50) by replacing 
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the elastic moduli and Poisson's ratios appearing on the right-hand 
side of the equation with their respective complex moduli and Poisson's 
ratios as follows: 


D 


12 


^*12 

- 2 r 3 E m R (l+i gEjn ) «[ J [F R (Tp+ i F 1 (T)) ] d T1 

o 

+ a/3)[(y6) 3 -i]v'(i+i Bvm )/n-v“ (1+1 8vm ) 2 ] } 


(171) 
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where F (t)) and F (T)) are the real and imaginary parts of the function 
F(T|) defined in equation (51) . Separating the real and imaginary parts 
on the right-hand side of equation (171), one readily obtains: 


°12 = 2 E ! F I Ol)]dTl +(1/3)[(u6) 3 -1]v R A R j- (172) 

0 o m 

dJ 2 = 2 E R r 3 | J [g E F R (n)+ F I (Tl)]dll + (l/3)[(U6) 3 -l]v R A 1 } (173) 

o m 

where 


A = (1+i g. ) (1+i g )/[l-v R (1+i g ) 2 ] 

E v m v 

mm in 

A R = Rea 1 ( A ) 

2 2 

“ ' (1-8. g )-(l+g p § )v (1+g )1 /A 
v E v Evm v ✓ A 

mm mm m 


A = Imaginary ( A ) 


2 2 

■ { <8 E + V )+<8 v ~ 8 E )V ! (1+8 v 3 } M 

mm mm; m 


> (174) 
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The loss tangent g^^ can now ca l cu l ate d from eqs. (L72-173) 
to yield: 


3 D 
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g D = { J Cgg F R (H) + F I (n)]drH-(l/3)[(U6) 3 -l]v^ A 1 } 

12 o m 


r 

l 




[F R (T0-g E F : (T 0 ]dW-(l/3)[(u6) 3 -l]v^ 

m 



(175) 


3.5 Transverse Flexural Loss Tangent 


Owing to antisymmetry in loading with respect to the T) axis 
as depicted in fig. 6, shear stresses t^ z again vanish as for the case 
of the transverse in-plane Young's modulus E ^ > Section 3.3. There- 
fore, the total strain energy and damping energies per cycle are given 

by the eqs. (165 and 166). The loss tangent g^>22 an< ^ t * le stora g e 
R 

stiffness are obtained by a procedure similar to that described 

in the first subsection of Section 3.3. 


3.6 Twisting Loss Tangent 


For a composite layer subjected to a pure twisting torque, all 

the stress components except t and t vanish. Thus, the strain 

xz yz 

energy and the damping energy per cycle in the fiber and matrix re- 
gions, respectively, are calculated from the following formulas: 
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(176) 


U 1 = IT (T 2 + T 2 ) dx dy / (2 G.) 

JJ xz yz 1 

A. 


U d ■ 2 " *0. u t 
1 


( i = f ,tn) (177) 


where the shear stress components and t are as given by eqs. 

(C-19 through C-22) for each typical composite element depicted in 

fig. 7. With the strain energies and damping energies calculated 

from eqs. (176 and 177), the loss tangent associated with the 

twisting stiffness D,, is obtained readily from eq . (151). 

b b 


3.7 Longitudinal Thickness-Shear Loss Tangent 

In view of eqs. (97) , the strain energies per cycle for the fiber 
and matrix regions are: 

1,1 = JJ (T xz +T yz )dX dy/(2 G i )+ JJ 

A. A. 

l l 


a 2 dx dy/(2 E.) 
z 1 

(i=f ,m) (178) 


where the strain components are as given by eqs. (97). According to 
eq. (147), the damping energies per cycle in the respective fiber and 
matrix regions are therefore equal to: 


I 

U d ' 2n g G. II (, L +T y 2 )dX d),/<2 G i ) 


+ 2 


" S E. II 


a dx dy /(2 E.) (i=f,m) 

z L 


(179) 
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The loss tangent associated with the longitudinal thickness- 

shear stiffness is then calculated according to the general procedure 
described in Section 3.1 with the help of eq . (151). 


3.8 Transverse Thickness-Shear Loss Tangent 

The strain energies per cycle for the fiber and matrix regions 
may be readily calculated from the results of elastic analysis in 
Section 2.8, eqs. (128, 129, and 135-137), as follows: 


- (9Q^/8G f )F 1 + { 3\ * Q^r4E m (u6) 6 ]}F 3 (180) 


= (3Q 2 /G ) {(u-1)/(5u6) + (3/8)F, } 

xml. ^ J 

+ {3Q 2 /[4E m <ii6) 3 ]} [(2/3)(u-l)(^i 2 + u+ 1) + F 4 ] (181) 


where the subscripts s and b refer to shear and bending, respectively; 
F^ and F£ are as defined in equations (136 and 137); and F^ and F^ are 
given by the following expressions: 



J {(u-|)/[l+(X , -l)(u6)‘ 3 (l-? 2 ) 3/2 ]} 2 (l-§ 2 ) 3/2 d^ 


(182) 
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(183) 


f 4 = J {<U-^/[i+(X'-l)(u6)' 3 (l-? 2 ) 3/2 ]} 2 

• [it6) 3 -(i-p 2 ) 3/2 ] 

In view of equations (181 and 182) , the damping energies per 
cycle in the respective fiber and matrix regions are equal to: 


= 2rtg G (9Q 2 /8G f )F 1 +2ng E {3\'Q 2 /[4E m <u6) 6 ]}f 3 (184) 

U d = 2TT8 G (3Q x /G m ) i (ii ' 1)/(5U 4 > + < 3/8 > F 2 } 
m 

+ 2t78 E { 3 ^/C 4E m (u6) 3 ]}[(2/3)(u-l)(^i 2 +u + l) +F 4 ] (185) 

m 

Finally, the loss tangent associated with the transverse thick- 
ness-shear stiffness is calculated on substitution of equations (180, 

2 

181, 184, and 185) into eq. (151). It is noted that , the square 

l 

of the shear force per unit z-length, cancels and thus the loss tangent 
is expressed in terms of the loss tangents, moduli ratios, and Poisson's 
ratios . 

This completes the damping analyses for the determinations of 
the loss tangents associated with all the stiffnesses pertinent to 
the characterization of the damping behaviors of a single layer of a 
monofilament composite. The results obtained herein are used in the 
next section. Section IV, for the construction of the design curves and 
comparisons of the numerical results with available analytical and 
experimental data. 
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SECTION IV 


NUMERICAL RESULTS 

In this section, numerical results as obtained from the foregoing 
elastic and damping analyses. Sections II and III, are presented. First, 
in Section 4.1, the results of the present elastic analyses are compared 
with some available analytical and experimental results obtained else- 
where as a verification of the present analyses. Secondly, in Section 4.2, 
property data deduction procedures for some of the constituent materials 
are briefly described. Then comparisons of the dynamic stiffnesses and 
their associated loss tangents are made with the limited experimental 
results which are available. Finally, in Section 4.3, the results of the 
present analyses are summarized in the form of a set of design curves for 
boron-epoxy composites, and in a tabulated form for boron-aluminum and 
glass-epoxy composites. This set of design curves is particularly useful 
in predicting the elastic and damping behavior of structural elements 
consisting of one or more layers of monofilament composites (reference 63). 

4.1 Comparison with Conventional Micromechanics Results 

This section is concerned with the numerical comparisons of the 
stiffnesses calculated from the present analyses with those obtained else- 
where in order to assess the accuracy of the present analyses. Two aspects 
of the comparisons to be considered here are: 
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1 . 


Comparisons of the in-plane and thickness-shear stiffnesses 
with those obtained from conventional micromechanics analysis 
where the composite layer consists of many small fibers randomly 
or regularly distributed throughout the entire cross section; 
see fig. 1(b) . 

2. Comparisons of the flexural and twisting stiffnesses with those 
deduced from the in-plane stiffnesses by the use of eq. (4). 

Major Young's modulus E^. ~ There have been many analytical in- 
vestigations to determine the major Young's modulus (refs. 1,39,47, and 
48) ranging from a simple mechanics-of -materials analysis to a more sophis- 
ticated elastic analysis. In all cases, however, it was demonstrated that 
the effect of difference in the Poisson's ratios of the constituent materials 
is negligibly small. In fact, Hill (ref. 47) showed by a variational method 
that the rule of mixtures, eq. (11), is the lower bound. Therefore, for all 
practical design purposes, eq. (11) is deemed sufficiently accurate. Further- 
more, experimentally, the published data for Narmco 3505 boron-epoxy com- 
posite (reference 64) gave values of to be 30.6 x 10^ psi in tension and 

34.0 x 10^ psi in compression; whereas, the value of 30.3 x 10^ psi was 
obtained from eq. (11) by using nominal Young's modulus values of 60 x 10^ psi 
for boron and 0.5 x 10^ psi for epoxy. 

Major in-plane Poisson's ratio - Many analytical investigations 

have shown that the major Poisson's ratio may be predicted by the rule of 
mixtures , eq. (38), with sufficient accuracy. In particular, Abolin'sh (ref. 
53) presented in detail the derivation of eq. (38) where the major Poisson's 
ratio of the fiber may be allowed to be transversely isotropic with the 
plane of isotropy normal to the fiber. However, in his derivation, no 
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allowance was made for the actual fiber cross-sectional shape. Within the 
context of mechanics-of-materia Is theory, Bert showed that eq . (38) holds 
for circular cross-section fibers (reference 65) . Halpin and Tsai (reference 
66) devised an interpolation method and applied it to the elastic numerical 
results in a graphical form obtained previously by Hermans (reference 67) to 
show that the numerical data can be empirically approximated by eq. (38). 

As a further verification of eq. (38), the numerical results for Vj^> as 
obtained from the second of eqs. (36), are shown in table I for various 
fiber volume fractions. It is observed that, for all practical purposes, 
the rule of mixtures, eq. (38), may be used for estimating the major Poisson's 
ratio v^ 2 * Unfortunately, comparison with the experimental result (ref. 64) 
showed that the theory predicted a much lower value of 0.30 against 0.36 
for Narmco 5505 boron-epoxy composite with a fiber volume fraction of 0.50. 
This is attributed mainly to the discrepancies in the nominal constituent- 
material data used for calculation. 

Minor Young's modulus E^ 2 ■ - The minor Young's modulus E 22 as 
calculated from the first of eqs. (36) and those obtained by other investi- 
gators are plotted against Young's modulus ratio E^/E^ for various fiber 
volume fractions for comparison as shown in figure 11. Foye has summarized 
the results of analytical investigations on the various estimates of the 
transverse properties of filamentary composites (ref. 40). It was concluded 
that, in general, the transverse stiffness is relatively insensitive to the 
types of models chosen. The specific analytical results for chosen 
here for comparison are those due to Tsai (reference 68) and Adams and Doner 
(ref. 56). In fig. 11, it is observed that the present theory predicts 
values for E 22 which are consistently lower than those of Tsai and of Adams 
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and Doner. However, it must be. remembered that in their analyses, a 
square-array composite model was used. Thus, the free-e.dge effect was 
completely eliminated. In the present analysis, a single-layer model 
was chosen in order to obtain a more realistic representation of the 
layer property on the micro-scale. Hulbert and Rybicki (reference 69) 
showed in their recent paper that the free-edge effect for some filamentary 
composites may range from 9.8% (for boron-epoxy) to 5.0% (for boron- 
aluminum) at a fiber volume fraction of 0.50. It is believed that this 
accounts for the lower estimates of the present analysis shown in fig. 11. Chen 
and Cheng (ref. 38) chose a hexagonal-array model and gave some numerical 
results for an E glass-epoxy composite and showed that their result from 
the elastic analysis is well within the previous results obtained by 
Hashin and Rosen (ref. 35) and Dow and Rosen (reference 70) . The present 
analysis also gave a result that is bounded by the results of the above 
two references (refs. 35 and 70). 

In-plane longitudinal shear stiffness G^. - Adams and Doner (ref. 

46) determined the in-plane shear modulus using a square-array model 
and compared their results obtained from an over-relaxation procedure with 
other analytical results (refs. 35, 38, and 71) and a limited number of 
experimental results. Excellent agreement was observed between their re- 
sults and the complex-var iab le elastic solution of Wilson and Goree (ref. 

71), whereas comparison with the analytical work of Chen and Cheng (ref. 

38) and that of Hashin and Rosen (ref. 35) showed some discrepancy due to 
the hexagonal-array model used by these investigators. A fair agreement 
is observed between Adam and Doner's results and that of experiment; 
theoretical values being consistenly lower by 5.37, for boron-epoxy, 6.7%> 
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for carbon-epoxy, and 13.5% for glass-epoxy composites. From an engineering 
design point of view, an explicit formula such as the one given by Hashin 
and Rosen ( ref. 35) would be highly desirable. For this reason, an 
empirical correction factor (tyf was introduced in Hashin and Rosen's formula 
in order to bring their results to coincide with that of Adams and Doner 
(ref. 46) . 

Longitudinal flexural stiffness D^. - In eq . (48), it is observed 
that the ratio represents a measure of flexural stiffness ef- 

ficiency which is less than unity for all realistic values of the parameters. 
In effect, this means that, unlike those composites containing many small 
fibers, the longitudinal flexural stiffness calculated from eq. (4) 
using the in-plane major Young's moudlus would be highly unconservative. 

A lower bound for the equivalent Young's modulus can be estimated readily 
from a "netting type" analysis in which the contribution of the matrix 
material to the composite flexural stiffness is completely neglected. This 
leads to the expression, originally derived by Margolin (reference 72). 


^• E ll ) /E ll-Wtti 


ing Analysis 


(3/4) /(a 5)' 


(186) 


In table II are shown the flexural stiffness efficiency, E^/E^, for 
various constituent-material combinations, aspect ratios 6, and volume 
fractions (reference 73) . It is noted that the effect of 6 on the flexural 
stiffness efficiency is much stronger than those of a and V^. 

Flexural stiffness efficiencies for various constituent-material 
combinations and fiber volume fractions are plotted as shown in figure 12 
for square typical elements, i.e., a square array of fibers. For example, 


75 


in the case of a boron-epoxy monofilament composite with a fiber volume 
fraction of 0.50, conventional theory, eq. (4), predicts a value for 
flexural stiffness that is twice as large as the actual flexural 
stiffness for a single layer. Also shown in the figure as dashed lines 
are lower-bound estimates for fiber volume fractions of 0.20 and 0.70. 

It is seen that the lower-bound estimate, eq . (186), increases in accuracy 
as the fiber becomes stiffer (large E^/E^ ratio), and as the fiber volume 
fraction increases. 

Poisson flexural stiffness - The values of Poisson flexural 

stiffness as calculated from the present analysis, eq. (50), are compared 
with those calculated from eq. (4) by the use of in-plane stiffnesses and 
Poisson ratios. Again, it is found that Poisson flexural stiffnesses cal- 
culated in a conventional fashion are much greater than those of present 
analysis as demonstrated in Table HI. For all of the composites compared, 

it is observed that for a low E../E ratio (<6) and low fiber volume fraction 

t m 

(<0.40), as calculated from eq. (4) is in fair agreement, the difference 
being less than 7%. However, for a high-stiffness-fiber composite such as 
boron-epoxy, the conventional formula, eq . (4), over-estimates by as much 
as 40% or more at a fiber volume fraction of 0.60. 

Transverse flexural stiffness D 22 • - In order to assess the trans- 
verse flexural stiffness efficiency in an analogous manner as that of 
longitudinal flexural stiffness, the ratio of the value of D 22 as obtained 
from eq. (59) and that obtained from eq. (4) is plotted against the E^/E m 
ratio for various fiber volume fractions as shown in figure 13. It is 
observed that the conventional estimates from eq. (4) are again on the 

unconservative side. However, variance in the transverse flexural rigidity 
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efficiency due to fiber volume fraction is relatively weak. This is 
perhaps due to the rather insignificant stiffening effect that the fiber 
has on the transverse in-plane and flexural stiffnesses. For a boron- 
epoxy composite with fiber volume fractions ranging from 0.4 to 0.6, 
the stiffening effect, as reflected by the value of ®22^m’ ran 8 es 
from 2 to 4, whereas in the case of the longitudinal modulus, the 

values of E../E ranges from 48 to 72! 

11 m 

Twisting stiffness D,,. - The torsional problem of a composite 

oo 

layer presented in Section 2.6 is particularized to a simple case of a 
single typical element and compared with the result obtained by Ely and 
Zienkiewicz, who solved the problem by the application of the relaxation 
method (ref. 74) . Comparisons of the values of the Prandtl torsion 
function at the mesh points showed that the present analysis is in good 
agreement with that of ref. 74; see figure 14. A series of exploratory 
computer runs indicated that the solution was relatively insensitive to 
the number of boundary points and the number of terms retained in the 
series solutions, eq. (80), for the values of parameters used in this 
example (G^/G^ = 10 » U " 2 , 6 ■ 1) . In the example shown, 31 equally 
spaced points on the vertical edge of the cross section and a 31-term 
series solution were used. 

In figure 15 is shown the values of the normalized Dirichlet 
torsion function v at equally-spaced mesh points of the one-quarter 
cross section of a composite layer containing three repeating typical 
elements. In this example, a total o I 50 points and 10-term and 20- 
term series solutions of the respective zeroth and first regions were 
used. Pertinent data are: G f /G = 0, u * 2, and 6=1. It is interesting 
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to note the similarity and symmetry in the V - value in the first and 
the second quarter-sections. It is also of interest to note the decrease 
in the V - value in the third quarter-section that is farthest from the 
origin. This means that for a layer containing many fibers, the 
torsional rigidity (not the layer twisting stiffness) can be approxi- 
mated as the number-of-f iber multiples of the torsional rigidity of a 
siqgle element. 

Figure 16 shows the twisting stiffness ratio for various 

values of fiber volume fractions and shear modulus ratio G^/G . 

f m 

Comparisons of the twisting stiffnesses for a square single 
element cross section calculated from eq. (84) and eq. (4) are tabulated 
as shown in Table IV. It is apparent that a twisting analysis is a must 
in predicting the layer-twisting stiffness. 

Longitudinal thickness-shear stiffness G,.^. - For a small-fiber 
composite or parallel laminates consisting of many layers, it is generally 
accepted that the longitudinal thickness-shear modulus G,-,. is equal to the 
in-plane longitudinal shear modulus G^ (references 75 and 76). However, 
in the case of a single-layer composite with only a single row of fibers, 
the longitudinal thickness-shear modulus G^ is expected to be greater 
than the in-plane shear modulus G^ due to the shear-stiffening effect of 
the fiber. Bert used an approximate Jourawski-type shear theory to pre- 
dict that the ratio G rr /G,, may vary from 2.86 for boron-aluminum to 37 
for boron-epoxy composites with a volume fraction of 0.482 (ref. 73). In 
the present more refined analysis, this ratio is found to vary from 1.5 
for boron-aluminum to 16.2 for boron-epoxy composites with a volume fraction 


of 0.5, 
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Numerical results for the longitudinal thickness-shear modulus 
are summarized as shown in figure 17 for various volume fractions and 
cons tituent-ma ter ia 1 combinations . 

Transverse thickness-shear modulus G... - The transverse th ick- 

44 

ness-shear modulus is plotted against volume fractions for various 

constituent-material combinations as shown in figure 18. Comparison with 
the analytical results obtained by Heaton (ref. 75) showed that fair 
agreement is observed for fiber volume fractions below 0.5, but at higher 
fiber volume fractions, 0.6 say, the present estimate for boron-epoxy 
composites gave a 14% higher value for G^. The discrepancy is attributed 
to the differences in the models chosen for the analyses; Heaton's model 
being that of a multi-fiber square array. 

4.2 Storage Moduli and Associated Loss Tangents. 

Prior to evaluations of the dynamic stiffness and damping be- 
havior of a layer of a monofilament composite characterized by nine 
stiffnesses and a Poisson ratio with their respective associated loss 
tangents, accurate dynamic constitutive properties of each constituent 
material must be known. A survey of the literature revealed that 
numerous experiments have been made to determine damping characteristics 
of glass (references 77-80), epoxy (ref. 9), and aluminum (references 
81-83). However, apparently no such data are available on boron. Because 
of different experimental techniques and perhaps slight differences in 
material compositions in the specimens, the damping properties obtained 
by these investigators do not always correlate with those obtained for 
the composite specimens. Therefore, it would be necessary to have avail- 
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able good data for the constituent-material properties for the micro- 
mechanics prediction of the composite macroscopic properties. However, 
sometimes this may be difficult to obtain experimentally, for example 
due to frailty of fibers in the filamentary composites or nonlinear 
effects (see Appendix D) . 

In view of this, a scheme for obtaining in-situ constituent 
material properties from tests on composites may be necessary. Some 
successful attempts at deducing such constituent-material properties 
were reported by Papirno and Slepetz (reference 84) and Bert (reference 
85) for static properties. Therefore, this deduction procedure was 
used in the present investigation properties to obtain the dynamic 
properties of some of the constituents from the test data on filamentary 
composites when necessary. 

To characterize the isotropic behavior of the material completely, 
two independent moduli and their associated loss tangents must be known. 
For the subsequent data deduction, it is assumed here that each material 
behaves elastically in dilatation . With this assumption and the knowledge 
of the Young's modulus and associated loss tangent, the second pertinent 
modulus (shear modulus or Poisson's ratio) and its associated loss 
tangent may be readily obtained. 

The dynamic Young's modulus and associated loss tangent for boron, 
epoxy, E glass, aluminum alloy 2024-T3, and aluminum alloy 6061 are 
summarized as shown in figures 19-23. Data for E glass (fig. 21) were 
obtained from ref. 79; for aluminum 2024-T3 (fig. 22), from ref. 81; 
for aluminum 6061 (fig. 23), from refs. 81 and 8; and curve B for epoxy 
(fig. 20) , from ref. 9. 
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Curve A for epoxy (fig. 20) was deduced from the experimental 
data on glass-epoxy filamentary composites by Mazza et al. (reference 
86); and the data for boron (fig. 19) were deduced from curve A and the 
experimental results on boron-epoxy filamentary composites of the 
same reference. 

A quick comparison showed that the prediction of the loss 
tangent associated with the transverse Young's modulus for glass-epoxy 
is in good agreement with the experimental results of ref. 9, both vary- 
ing from 1% to 3 °L for the frequency range between 20 and 2000 Hz. 

The same comparison for boron-epoxy filamentary composites also 
indicated excellent agreement on the prediction of the loss tangent 
associated with the transverse Young's modulus; in both cases, ref. 86 
and the present analysis, the value of loss tangent varies between 1.6 
and 2.07o in the frequency range 20-400 Hz. 

Pertinent data for the complete characterization of the con- 
stituent properties are summarized in Tables V-IX. 

4.3 Design Curves and Tables - Storage Moduli and Associated Loss 
Tangents Versus Frequency for Various Fiber Volume Fractions 

The results of the present analyses are summarized in the form 
of design curves for a boron-epoxy filamentary composite as shown in 
figs. 24-33 for the frequency range 50-2000 Hz. and the fiber volume 
fractions V f = 0.4, 0.5, and 0.6. In application, the curve A series 
should be used for better predictions of the composite stiffnesses 
and associated loss tangents; the Curve B series is included for com- 


parison purposes only. 
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In Tables X-XX1X are listed all the pertinent dynamic stiffnesses 
and associated loss tangents for the characterization of the layer 
properties of boron-aluminum and E glass-epoxy composites for the 
same ranges of frequency and fiber volume fractions. 


V. CONCLUSIONS 


Elastic and damping analyses resulting in determination of 
all pertinent stiffnesses and associated loss tangents for the 
characterization of the elastic and damping behavior of a mono- 
filament composite were carried out. 

The numerical results obtained for the stiffnesses and 
associated loss tangents compared favorably with some existing 
analytical and experimental results for some typical filamentary 
composites, such as, boron-epoxy, boron-aluminum, E glass-epoxy. 

The results of the flexural and twisting stiffness analyses 
showed that these properties cannot be deduced accurately from the 
in-plane-properties, and, thereby, the necessity for such analyses 
for a monofilament composite layer. 

The assumption of Kimba 1 1-Lovel 1 type damping was shown to 
be equivalent to the elastic-viscoelastic correspondence principle 
for the case of the in-plane longitudinal stiffness. 

The results of this investigation were summarized in a set 
of design curves for a boron-epoxy composite, and in a set of design 
data tables for boron-aluminum and E glass-epoxy composites. The 
former were applied to the problem of predicting resonant frequencies 
nodal patterns, and damping ratios for laminated boron-epoxy plates 
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(references 64 and 87 ) and achieved excellent agreement with available 
experimental results. for six different plates (ref. 14 ) . 

Recommendations for future investigation . - For future re- 
searches, the following topics are suggested: 

1. Experimental characterization of the complete set of dynamic 
stiffness and damping properties of various filamentary 
composite materials of technological importance. 

2. Analytical investigation of the viscoelastic behavior of 
composite materials with thermo-mechanical coupling. 

3. A unified viscoelastic analysis for filamentary composite. 
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APPENDIX A 


BOUNDARY-POINT LEAST-SQUARE METHOD 

Numerous approximate methods are available for the solution of 
boundary-value problems. To name a few, there are the Rayleigh-Ritz , 
Galerkin, Kantorovich, finite element, relaxation, and collocation 
methods. Of these methods, especially with the availability of modern 
digital computers, probably the boundary- point least-square version 
(ref. 5) of the collocation method is the most efficient for solution of 
mixed boundary-value problems, such as those investigated in Section 
II. The main useful features of the method are: 

1. It may be applied to mixed boundary-value problems with re- 
lative ease. 

2. The assumed solution may be made to satisfy the boundary conditions 
at a set of sufficiently dense points in the sense of minimiz- 
ing the square error; hence the solution may be made independent 

of the number of boundary points chosen. 

In general then, the solution is reduced to the satisfaction of a set 
of overdetermined algebraic simultaneous equations 

AX = B (A- 1) 

where m x n coef f ic ient matrix (m>n), X -- n-dimension column vector 

of unknown coefficients, B = n-dimension column vector of prescribed 
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boundary values, m = number of chosen boundary points, and n = number 
of unknown coefficients. 

2 

Then the mean-square-error matrix (E ) for equation (A-l) is given 
as 

0 ****** nr ******* 

E = (AX-B) (AX-B) (A -2) 

where the superscript T denotes the transpose of the matrix quantity 
which it follows. 

On minimizing equation (A-2) with respect to X, one obtains 


(X T A)X = (a t b) 


(A-3) 


or 


A X = B 


(A-4) 


where A is aa n x n matrix given by 


~ * 

A 


~T~ 

A A 


and B is 


an n-dimension column vector given by 


.-Jf 

B 


-S t b 


(A -5) 


(A-6) 


Finally, the solution-coefficient vector X that minimizes the square 
error is obtained readily from equation (A-4) by using a number of standard 
computer scientific subroutines. However, care must be taken in the choice 
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of appropriate forms of the stress functions and in the choice of 
boundary points to avoid any boundary points at which the prescribed 
boundary values are identically satisfied. 



APPENDIX B 


MODELS AND MEASURES OF MATERIAL DAMPING 
Bl. Mathematical Models for Material Damping 

Internal friction or damping in materials can be caused by a variety 
of combinations of fundamental physical mechanisms, depending upon the 
specific material. For metals, these mechanisms include thermoelasticity 
on both the micro and macro scales, grain boundary viscosity, point- 
defect relaxations, eddy-current effects, stress- induced ordering, inter- 
stitial impurities, and electronic effects (ref. 7). According to 
Lazan (ref. 7 ), little is known about the physical micromechanisms 
operative in most nonmetallic materials. However, for one important class 
of these, namely polymers and elastomers, considerable phenomenological 
data have been obtained. Due to the long-range molecular order associated 
with their giant molecules, polymers exhibit rheological behavior inter- 
mediate between that of a crystalline solid and a simple liquid. Of 
particular importance is the marked dependence of both stiffness and damp- 
ing on frequency and temperature. 

The purpose of developing a mathematical model for the rheological 
behavior of a solid is to permit realistic results to be obtained from 
mathematical analyses of complicated structures under various conditions, 
such as sinusoidal, random, and transient loading. According to reference 
7 , as early as 1784 Coulomb recognized that the mechanisms of damping 
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operative at low stresses may be different than those at high stresses. 

Even today, after more than 2500 publications on damping have appeared, 
major emphasis is placed on linear models of damping for several reasons: 
they have sufficient accuracy for the low-stress regime and linear analyses 
are computationally more economical than nonlinear ones. 

The simplest mathematical models of rheological systems are single- 
parameter models: (1) an idealized spring, which exhibits a restoring 

force linearly proportional to displacement and thus displays no damp- 
ing whatsoever, and (2) an idealized dashpot, which produces a force 
linearly proportional to velocity. Obviously, neither of these models 
are very appropriate to represent the behavior of most real materials. 

The next most complicated models of rheological systems are the 
two-parameter models: (1) the Maxwell model, which consists of a spring 

and dashpot in mechanical series, as shown schematically in figure B-l(a); 
and (2) the Kelvin-Voigt model, which is comprised of a spring in 
parallel with a dashpot, as shown in figure B-l(b). The Maxwell model 
is a fair approximation to the behavior of a viscoelastic liquid. How- 
ever, as a model for a viscoelastic solid, it has several very serious 
drawbacks (ref. 7) ; there are no means to provide for internal 
stress and for afterworking. The Kelvin-Voigt model overcomes these 
deficiencies and is a first approximation to the behavior of a viscoelastic 
solid. However, it has the following disadvantages (ref. 7): there 
is no elastic response during application or release of loading, the 
creep rate approaches zero for long durations of loading, and there is 
no permanent set irrespective of the loading history. 
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It should be mentioned that the Kelvin-Voigt model is the simplest 
one which permits representation as a complex quantity when subjected to 
sinusoidal motion. To show this, one can write the following differential 
equation governing the motion of a single-degree-of-freedom system con- 
sisting of a sinusoidally excited mass attached to a Kelvin-Voigt element: 

mu + cu + k u = Fe (B-l) 

where m s mass, c = viscous damping coefficeint, k = spring rate, F h 
exciting force amplitude, iu = circular frequency of exciting force, 
u = displacement, i ^/~T ~ , t s time, and a dot denotes a derivative 
with respect to time. Assuming that sufficient time has passed for the 
transients to die out, we can write the steady state solution of equation 
(B-l) as follows: 

u = u e i ^ U>t 1" [(k-muj 2 ) + iu)c] 1 e iU)t (B-2) 

where cp = phase angle between the response (u) and the exciting force, 
and u & displacement amplitude. 

It is noted that the terms containing Kelvin-Voigt coefficients k and 
iu)c can be grouped into a single Kelvin-Voigt complex stiffness (k) as 
follows : 


It can be shown that they do die out for a mechanical system, for which 
c, m, and k are all positive quantities. 
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k = k + iujc 


(B-3) 


Alternatively, the Kelvin-Voigt element coefficient could have 
been grouped into a single complex damping coefficient (c ) as follows: 

c = c - (i k/oj) (B-4) 

The force (F^) in a dashpot is given by 

„ ~ i(u)t - cp) C \ 

Fj = c u = tujcue ^ (B-5) 

The energy dissipated per cycle (U^) in a dashpot undergoing 

sinusoidal motion can be calculated as follows: 

2rr/u) ~ 

F, u dt = TTCivu (B-6) 

0 d 

In 1927, Kimball and Lovell (ref. 22) observed that many 

engineering materials exhibited energy losses which were contradictory 

to equation (B-6). Specifically, they found that was proportional 
~2 

to u but independent of oj, i. e. 

U d = Cju 2 (B-7) 

~2 

Later work by Wegel and Walther (reference J8) also showed that « u 

but Cl was a weak function of w- Still later, for stresses below the 
a 

fatigue limit of the material, Lazan (references 7,89 ) showed that « 
and Cj was practically independent of uj. 


= i 


F, du = 


Equating the energy dissipated per cycle in the Kimbal 1-Love 11 
material, equation (B-7), to that in the dashpot, equation (B-6) , one 
obtains 

C* = new (B-8) 

d 

To utilize the Kimbal 1-Love 11 observation in practical structural 
dynamic analyses, it was desirable to have an expression for the as- 
sociated damping force, i.e. an equation analogous to equation (B-5) , 

I 

which governs a dashpot. According to Bishop (reference 90) this was 
first done by Collar. Combining equations (B-5, B-8), one obtains the 
following "frequency-dependent damping" relation: 

F d = (b/o>) u (B-9) 

where b is a constant given by 

b = C^/tt (B- 10) 

The kind of damping represented by equation (B-9) has been called fre- 
quency-dependent damping, since the usual dashpot coefficient c is now 
replaced by the quantity b/w. 

For the damper represented by Eq. (B-9), it is convenient to replace 
k in the undamped system by the Kimball-Lovell complex stiffness 
(reference 91): 

k' = k + ib (B— 11) 

where it is now noted that both k and b are assumed to be independent of 
frequency. This representation has been used extensively in aircraft 
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structural dynamic and flutter analyses (references 92,93). 

An alternative to equation (B-ll), is to use the concept of a complex 
damping coefficient, originated by Myklestad (reference 94), in which the 
spring constant is replaced by 

k' ' - Cl e Lm (B-12) 

where and m are constants. 

Although Kimball and Lovell's original work was based on data 
obtained during steady-state sinusoidal motion, the models represented 
by equations (B-9, B-ll, B-12) have been applied to free vibrations as 
well (refs. 90, 91 ,95). To alleviate the difficulty of the interpreting of 
u) in equation (B-9) for the case of free vibration or for multiple-fre- 
quency forced vibration, Reid (reference 96) suggested the following form: 

“ b | u/u j u (B-13) 

where the bars denote that the absolute value of the quantity between 
them is to be used. 

In summary, four different versions of material damping have been 
discussed. In terms of the differential equation for free vibration, 
they may be written as follows: 


mu 

+ (b/uj) u + ku 

= 0 

(B- 14a) 

mu 

+ (k + ib) u 

= 0 

(B- 14b) 

mli 

„ im 
+ C^e u 

= 0 

(B- 14c) 

mu 

+ b I u/u j u + 

ku = 0 

(B- 14d) 
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Tne models represented by equations (B-14) have been criticized 
by many investigators for various reasons: 

1. In equation (B-14a), the interpretation of uo is dubious (refer- 
ences 96,98). Milne (reference 97) suggested that it is associated with the 
imaginary part of the pair of complex roots, but went on to state that 
there is no clear justification for this interpretation. 

2. Equations (B-14b) and (B-14c) are equations with complex coef- 
ficients and the meaning of this is not clear (refs. 95, 9/). Further- 
more, although they have complex exponential solutions, neither the real 
nor imaginary parts alone are solutions (refs. 97, 98). 

3. Equation (B-14d) is a nonlinear differential equation, due to 
the behavior of | u/u| (ref. 97). 

4. None of the models represented by equations (B-14) can be sim- 
ulated, even on a conceptual basis, on an analog computer (refs. 

96, 98 ) . 

5. The model represented by equation (B-14a) does not meet the 
Wiener-Paley condition (references 99,100) of causality for physically 
realizable systems, as was shown in references 98,102-103. 

6. Equation (B-14b) fails to give the proper relations of a 
damping force that is proportional to displacement and in phase with 
velocity. This can be alleviated by use of a more unwieldy expression 
proposed in references 104,105. 

Inc idental ly ref erence 106 claimed that none of the equation 
(B-14) models result in energy dissipation that is independent of 
u) (the reason for abandoning the Kelvin-Voigt model) . However, an error 
in this reasoning, which makes it invalid, was pointed out in reference 107. 
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It can be shown that the sinusoidally forced vibration version of 


free vibration equations (B-14a,b,c) are exactly equivalent. Probably 
the most popular form for writing this is as follows: 


mii + (k + ib)u 


-i, iu>t 
Fe 


(B- 15) 


The major criticism of the use of equation (B-15) is that it results 
in a conventionally defined magnification factor (see Section B2), correspond- 
ing to zero frequency, which is not equal to unity for g ^ 0 (ref. 90) . 

This difficulty can be eliminated by redefining the magnification factor as 
follows : 

MF' = u/u^ (B- 16) 


where 

u' = *F/k' 
st 

and k' is the stiffness coefficient associated with the total in-plane force 
amplitude, i.e. 

k' s (k 2 + b 2 )^ 


Thus , 

MF' = (k 2 + b 2 A'[k-nxu 2 ) 2 + b 2 ]^ 


2 ,^ /r . 2 . 2,2 2 ^\ 
(1 + g ) /[ (1 -oj /a> n ) + g ] 


The major advantages of using the equation (B-15) model for steady- 
state vibrational and flutter analyses are as follows: 

1. The analysis is considerably simplified in comparison with that 
for viscous damping, as pointed out in reference 108. 

2. The analysis is linear. 

3. The complex-modulus concept has long been in common use, especially 


95 



in connection with measurement of the damping properties of elastomers, 
polymers, etc. (reference 107). 

4. The complex-stiffness approach has been used very successfully 
in flutter analysis for approximately three decades (refs. 92,93 ,107 ,109 ). 

In view of these considerations, the complex-stiffness approach is 
used in the main portion of the present report. However, for completeness 
and comparison, a number of other more complicated approaches are also dis- 
cussed in this section. 

To overcome some of the deficiencies of the Kelvin-Voigt and Max- 
well models, they were combined in various ways to obtain the three-para- 
meter models shown in figure B-l(c) and (d). It can be shown (reference 
110) that by properly selecting the values of the coefficients, the two 
models can be used interchangably . In other words, the model repre- 
sentation is not unique. In fact, both representations have been called 
the "standard linear solid", "standard linear material," "simple an- 
elastic model", or "standard model of a viscoelastic body" (references 
111-114)* Unfortunately, few materials have been characterized by the use 
of the standard linear solid. Perhaps the most extensive characterization 
has been carried out for flexural vibration of 2024-T4 aluminum alloy: 
by Granick and Stern (reference 82,115) in air and by Gustafson et al. 

(ref. 81) in vacuo . 

Obviously one can go on from a three-parameter model to a four- 
parameter one, etc. Continuing this process, one obtains two well- 
known models (ref. Ill): the Kelvin chain, shown in figure B-2(a), 

and the generalized Maxwell model , shown in figure B-2(b). The be- 
havior of such systems can be written in either differential or integral 
form, as discussed in detail in reference 111. However, the complexity of 
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such a model has limited its use both experimentally (in characterizing 
engineering materials) and analytically (in performing dynamic analyses). 
Thus, for engineering applications, there is a continuing search for 
models which are simple, yet represent the behavior of real materials 
with sufficient accuracy. 

A relatively simple model, which is somewhat similar to the generalized 
Maxwell one, is due to Biot (reference 116); see figure B-2(c). The number 
of simple Voigt elements in this model is allowed to increase without bound, 
so that the damping force is represented by the following expression: 


pt 


®1 


u t 

where t is a dummy variable; e , 
integral 



[ -e (t-T )] (du/dr) dr 
parameters, and Ei(u) 


(B-17) 

is the exponential 


Ei (u) 5 (e' ? /§) d§ (B-18) 

Caughey (ref. 98) made a detailed study of Biot’s model for both 
free and forced vibration. He showed that for u)/e > 10, the energy dis- 
sipated in Biot's model is within 4% of being independent of frequency (to) 
and thus essentially depends only upon the square of the displacement 
amplitude. Milne (ref. 97) studies the Biot model by means of the con- 
volution integral. 

Neubert (reference 117) introduced a variant of the Biot model, in 
which a different distribution function for the local stiffness-damping 
ratio is used. Milne (ref. 97) also studied this model via the convolution 
integral. Milne introduced a new synthesized model of his own. However, 
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the main conclusion of his investigation was that, for practical engineer- 
ing purposes, the choice of a rheological model is not critical, provided 
that the range of constant damping is kept the same and the damping is 
sma 1 1 . 

The linear hereditary theory of material damping was originated in 
1876 by Boltzmann (reference 118) . In this theory, the energy loss is 
attributed to the elastic delay by which the deformation lags behind 
tne applied force. It is called a hereditary theory because the in- 
stantaneous deformation depends upon all of the stresses applied to the 
body previously as well as upon the stress at that instant. The damping 
force is given by 

F d = f §(t,-r) u(r) d T (B- 19) 


where t = actual time, T = an instant of time (dummy variable), u = 
displacement, and $ = hereditary kernel or memory function. Often the 
hereditary kernel depends upon the difference (t-T) only; then equation 
(B-19) becomes: 


I 


t 

$ ( t - t) u(t) dr 

-00 


(B-20) 


It is most advantageous to determine the hereditary kernel from 
experimental data. It has been found to be a monotonically decreasing 
function which can be represented mathematically as follows: 

n 

§(t) = £ A.e" Qf i t (B-21) 

i = 1 

Such a system is said to have a heredity of degree n and and 
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are called heredity constants. Hobbs (ref. 112) has shown that a first 
degree hereditary model is equivalent to a simple anelastic model 
(standard linear solid) . 

Volterra (reference 119) has presented a method of determining the 
hereditary kernel from material response to either suddenly applied 
loading or to sinusoidal excitation. He also presented plots of mag- 
nification factor and phase shift versus frequency for a wide range of 
parameters of a single-degree-of- freedom system with first-degree 
hereditary damping. For certain combinations of the parameters, the 
magnification and phase characteristics are quite different from that 
of a viscously damped system; yet for certain other combinations, the 
characteristics are quite similar. 

Numerous nonlinear mathematical models have been proposed to per- 
mit better correspondence between theory and experiment. However, all 
of them, by their very nonlinear nature, are more complicated to apply 
in structural dynamic analyses of practical engineering systems. Perhaps 
the simplest one is damping energy per unit volume and per cycle pro- 
portional to stress to a power (ref. 7 ): 

U dv - c~a m (B- 22) 

"V 

where and m are material constants and a is the cyclic stress 
\amplitude. It is noted that m = 2 corresponds to the Kimball-Lovell re- 
lation (linear theory) discussed previously. It has been found that m 
varies from 1.8 to 6.0, generally being near 2.0 at low stress amplitudes 

Tfte hysteresis loop associated with the complex-modulus model, 
equation (B-ll), is elliptic in shape; see figure B-3(a). However, 


the hysteresis loops determined experimentally have sharp corners at both 
ends and are nearly linear in the low-stress, low-strain region, as shown 
in figure B-3(b). According to ref. 7 (page 97), in 1938 Davidenkov 
proposed a nonlinear mathematical model which results in a similar shape 
of hysteresis loop. Further work on this class of model has been carried 
out by Pisarenko (reference 120) . 

Rosenblueth and Herrera (reference 121) proposed a nonlinear model 
which eliminates the causality difficulties of the complex-modulus model. 
Chang and Bieber (reference 122) introduced a nonlinear hysteresis model 
in which there is no hysteresis damping unless the displacement amplitude 
exceeds a certain threshold value. 

B2 . Measures of Material Damping 


In this section, various definitions of damping are discussed in 
the context of a complex-modulus material. 

Energy Dissipation Under Steady-State Sinusoidal Vibration. - Energy 
dissipated per cycle (U^) in the form of internal frictional heating is 
one measure of damping. However, this quantity depends upon the size, shape, 
and dynamic stress distribution (which in turn, depends upon the particular 
mode of vibration). In view of the above difficulties, the specific damping 
energy (U^ ) is usually considered to be a more basic property of the 
material , rather than the structure. The U^ y is defined as the damping energy 
per cycle and per unit volume, assuming a uniform dynamic stress distribution 
throughout the volume considered. Thus, the total damping energy is 


U 


d 



dv 


dV 


( B- 23) 
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where V = volume. The usual units of U, are in-lb/cvcle and of U, , 

a dv 

3 

in-lb/in -cycle. 

Resonant Magnification Factor Under Steady-State Sinusoidal Ex- 
citation . - The resonant magnification factor (RMF) is defined as the 
dimensionless ratio of the response at resonance to the excitation, 
where both the response and the excitation must be specified in the 
same units (see figure B-4) . These quantities may be in units of dis- 
placement, velocity, acceleration, or strain. Unfortunately, the re- 
sonant magnification factor is dependent upon the structural system con- 
figuration as well as upon the damping property of the material, so that 
it is considered to be a system characteristic, rather than a basic material 
property. 

Incidentally RMF is not applicable to nonlinear systems, since then 
it depends upon the excitation level. (In a linear system, RMF is independ- 
ent of the level of excitation.) 

Bandwidth of Half-Power Points Under Steady-State Sinusoidal Excitation . - 
The separation (uj^-uj^) between the frequencies associated with the half- 
power points increases with an increase in damping (see figure B-5) and thus 
can be used as a measure of damping. A more meaningful definition is to 
normalize by the associated resonant frequency , i.e. use the dimension - 

less bandwidth given by 

(u) 2 - UJ 1 )/ UJ n 

This measure is the basis for determination of damping by the original 
Kennedy-Pancu method (reference 123) . 

The quality factor Q is defined as the reciprocal of the dimensionless 
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bandwidth. 


Derivative of Phase Angle with Respect to Frequency, at Resonance . - 
This measure is the basis for determination of damping by an improvement 
(reference 124) of the Kennedy-Pancu method. 

Loss Tangent Under Steady-State Sinusoidal Excitation . - Applying 
the complex-stiffness (see Section Bl) to the material property represent- 
ing stiffness, namely the Young's modulus (E) , we obtain 

E - E R (1 + ig) (B-24) 

ic 

Thus, the loss tangent is defined as follows: 

g = E X /E R (B-25) 

j ** g *** 

where E = loss modulus and E is called the storage modulus . 

Referring to figure B-6, we obtain the following relation: 

g = tan y = E*/E R (B-26) 

or 

y = arctan E*/E R (B-27) 


*Also known as the "loss coefficient", "loss factor", or "damping 
factor". 

**Also sometimes called the "dissipation modulus." 

***Also sometimes called "elastic modulus" or "real modulus." 
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The quantity y is called the loss angle and equation (B-26) explains 


the origin of the term loss tangent for g. 

Cyclic Decay of Free Vibrations . - In a linear system, free vibrations 
decay exponentially (see figure B-7). The larger the damping, the faster 
is the decay. Thus, the logarithmic decrement 6 is defined as follows: 

6 s In (a /a ) (B-28) 


where £n s natural logarithm. 

For a power-law material, eq. (B-22), to have 6 independent of 
amplitude, m must be 2 (see Section B3 for proof). When m=2 , the follow- 
ing means of calculating 6 is more practical, especially when damping 
is small: 


6 = (1/n) -tn (a / a ) (B-29) 

1 l-Hl 

where n is any arbitrary integer. 

Temporal Decay of Free Vibrations . - Another measure associated with 
the decay of free vibrations is the temporal decay constant v^, defined as 
follows (reference 7 , foldout opposite p. 35): 

V t ^ ( t 2 ~ t p 1 C u(t 1 )/u(t 2 )] (B-30) 

where t^ and are two different values of time, arbitrary except that 
t 2 > t^ and (t 2 - t^) must be an integer number of periods of damped 
vibration. The quantities u(t^) and u(t 2 ) are the respective displacements 
at times t^ and t 2> It is seen that the usual units for v^_ are sec 
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An alternate way of specifying the temporal decay is the decay 


rate y^, which has units of decibels per second (db/sec) and is defined 
as follows: 

Y t = 20 (t 2 - t L ) _1 log [u( tl )/u(t 2 )] (B-31) 

or 

-1 V t (t 2 " t l ) 

Y t = 2 °(t 2 " t^) log e = (20 log e)v t 8.68 v fc 

where log s logarithm to the base 10. 

Spatial Attenuation of a Plane Wave in a Slender Bar . - If 
one propagates a plane, harmonic wave in a slender bar made of a homo- 
geneous, linear material, the wave exhibits exponential decay with axial 
position (x), analogous to the temporal decay of free vibration (which 
may be considered as a standing wave). Analogous to the logarithmic de- 
crement (see above), we have the logarithmic attenuation 6 , defined as 

s 

follows in terms of amplitudes a(x^) and a(x 2 ) at stations x^ and x 2 : 

6 g s -tn [a(x 1 )/a(x 2 )] (B-32) 

Analogous to the temporal decay constant, we define the spatial 

attenuation constant v as follows: 

s 

V g = (x 2 - Xj^) 1 tn r a ( x ^)/ a (* 2 )] (B-33) 

The unit of v is in 
s 

Finally, analogous to the decay rate, we define the spatial 

attenuation rate v as follows: 

Y s 
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(B-34) 


Y s s 20 (x 2 - 1 log [a (x^/a (x^)] % 8.68 v g . 

The unit of Y g is db per unit length. 

It should be cautioned that the determination of damping by wave 
attenuation requires that the cross section of the bar be compact (pre- 
ferably round or square) and that the largest cross-sectional dimension 
be very small compared to the wave length of the traveling wave. Further- 
more, the two stations along the bar, at which the wave measurements are 
made, must be: (1) sufficiently far from the point of impact that initial 
transients have died out and (2) sufficiently far from each other so that 
a measurable change in amplitude occurs. Kolsky discussed these points 
in detail in reference 125. 

B3. Inter-Relationships Among Various Measures of Damping for 

Homogeneous Materials 

Here we consider a single degree-of- freedom system consisting of a 
mass supported on a massless spring made of a homogeneous Kimball -Love 11 
or complex-modulus material. Then the governing differential equation 
for simple harmonic excitation is equation (B-15), which has the following 
steady-state solution: 

~ i(u)t-o) -jc\ 

u = u e Y (B-35) 

where 

u =s F [ (k-muj 2 ) 2 +b 2 ]" 1/2 (B-36) 

cp = arctan [b/ (k-moj 2 )] (B-37) 

The magnification factor (MF) is defined as follows: 
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MF s u/u 


(B-38) 


st 

where -u is the so-called static displacement, defined as follows: 

u =J/k (B-39) 

st 

Combining equations (B-36, B-38, and B-39), one obtains the 
following relation: 

MF = {r 1 " (u>/up n )2]^+ (b/mu) 2 ) 2 J- *^ 2 (3-40) 

where u)^ is the natural frequency (i.e. the resonant frequency of the 
undamped system), given by: 

tt) n = (k/m) 1/2 (B-41) 

The critical material damping coefficient (b c ) is the minimum value 
of b for which the free vibratory motion is non-oscillatory. It is 
obtained from the following relationship: 

(b /2mu) ) 2 = k/m = uj 2 (B-42) 

c n n 

or 

b = 2mu> 2 (B-42) 

c n ' 

The damping ratio £ is defined as follows: 

C s b/b c (B-43) 

Combining equations (B-40, B-42, and B-43), one arrives at the 
following general, dimensionless relationship: 
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MF = {[l-(u>/a> n ) 2 ] 2 + 4£ 2 }" 1/2 


(B-44) 


It Is seen easily that MF is a function of only two parameters: the fre- 

quency ratio (u>/u> n ) and the damping ratio (Q ) . It should be mentioned 
that equation (B-44) is somewhat different than the known relationship 
for a s ingle-degree-of-freedom Kelvin-Voigt viscously damped system, 
even though the damping ratio Q is defined analogously ( = c/c^, 

where c = 2mjj ). 
c n 

The resonant magnification factor (RMF) is defined as the value 
of MF at resonance, which is taken here to occur when cp = 90°. In view 
of equation (B-37), this implies that the resonant frequency is 
Thus, at resonance, equation (B-44) reduces to the following simple 
expression: 


RMF = (20 


(B-45) 


Thus, it is seen that the same simple relationship between the damping 
ratio and the resonant magnification factor holds as does for a Kelvin- 
Voigt system. Although the resonant amplitude can be measured accurately, 
there are problems in measuring the static deflection in complicated 
structures (due to signa 1-to-noise ratio limitations of instrumentation). 
For this reason, in structural dynamics, RMF is seldom used as a measure 
of damping. 

The damping energy, i.e. the energy dissipated per cycle, can be 
computed as follows: 

2rr/u) 

o p 

U d = j F d du = I F d a dt (B-46) 

where F^ is the damping force, given by 
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F d = (b/u))u 


(B-47) 


Substituting equations (B-35) and (B-47) into equation (B-46) and 
integrating, we obtain the following result: 

U, = nb u 2 (B-48) 

d 

To convert equation (B-48) to the form of equation (B-22), with 
m=2 , it is necessary to divide by the volume and to relate u to < 7 . 

Now the force amplitude is 

F = /£ = ku (B-49) 

where A h cross-sectional area. 

Assuming the spring is a massless and uniform bar, we have 

k = AE/L (B-50) 

where L = effective length of spring and E = Young's modulus of the bar. 
Combining equations (B-48, B-49, B-50), one obtains the following result: 

U dv = U d /V = Cd bIj / 4) (o/E) 2 (B-51) 

Comparing the form of equations (B-22) and (B-51), we see that 
for a Kimball -Love 11 material, m = 2 and 

C d = nbL/AE 2 (B-52 ) 

It should be noted that Lazan (reference 7 , chap. II) has shown 
that these results are independent of the stress distribution and the 
specimen geometry, provided that m = 2. 
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In view of equations (B-42 , B-43, and B-44) , equation (B-52) can 
be expressed in the following more useful form: 

(B-53) 

According to equation (B-48), the damping energy (and thus, 
excitation power) is proportional to the square of the displacement 
amplitude u q (and thus, the MF) . Then the power is one half of what it 
is at resonance when MF = RMF//2. In view of equation (B-45), the half- 
power magnification factor (HPMF) is related to the damping ratio as 
follows : 

HPMF = RMF//2 = (2/2 £)' 1 (B-54) 

To determine the two sideband frequencies (id = 10 ^ and cu^) cor- 
responding to half power, we combine equations (B-44) and (B-54) to 

2 

obtain the following quadratic expression in id : 

(2 72 C ) 2 = [1 - (w/w n ) 2 ] 2 + 4£ 2 
which has the following positive roots: 

u>l 2 = u, n (1 - 2 (B_55) 

It is interesting to note that equation (B-55) is identical to the analogous 

2 

expression for the lightly damped Kelvin-Voigt system, i.e. ( c/c ) «1. 

For small damping (Q «%) , 

U>1 2 (1 ± C) 
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Thus, we get the following relation between the dimensionless bandwidth 
and the damping ratio Q: 

(B-56 ) 

In the original Kennedy -Pancu method (reference 123), this relationship is 
used to determine the damping ratio from the geometry of an experimentally 
determined Argand plot (a polar-coordinate plot of response amplitude 
versus phase angle). 

In terms of the quality factor Q, equation (B-56) may be rewritten 
as follows: 

(B-57) 

It is interesting to note that in order for the half-power frequencies 
to be real, the damping ratio Q must be no greater than /2/2 s a 0.707. How- 
ever, this is not a severe limitation for actual structural materials. 

The potential energy (strain energy) stored in our s ingle-degree-of- 
freedom mathematical model is given by: 

/■> 

U = o F g du (B-58) 

where F is the spring force, given by 
s 

F = ku (B-59) 

s 

Thus , 

u 

U = u du = ku^/2 (B-60) 

Combining equations (B-48) and (B-60), we obtain 




U d /U = 2nb/k 


(B-61) 


no 




Using equations (B-41, B-42 , B-43), we can rewrite equation (B-61) 
in the following more useful form: 

U d /U = 4 ttC (B-62) 

or 

(B-63) 

Hie loss tangent is defined as follows for the system considered 

here : 

g = b/k (B-64) 

Inserting equations (B-41, B-42, B-43) into equation (B-64) we 
obtain the following useful relationship: 

(B-65) 

The following useful relationship, which is sometimes taken as a 
fundamental definition of damping (ref. 104), is obtained here by 
combining equations (B-63) and (B-65): 





At resonance (subscript R) , defined by u> = u> n ^/k/m, equation (B-67) 
becomes : 

(dcp/duu)^ = 2mu) n /b (B-68) 

Combining equations (B-42, B-43, and B-68), we obtain the following 
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very useful result: 


c = [u> (dcp/du)) ] 1 > (B-69) 

It is interesting to note that although the MF vs. u) relation, 
equation (B-44), is different than the corresponding relation for a Kelvin- 
Voigt viscously damped system, equation (B-69) is identical to its cor- 
responding relationship for a Kelvin-Voigt system. Equation (B-69) is used 
to determine the damping ratio from an experimentally determined Argand 
diagram in the improved Kennedy-Pancu method (ref. 124) . 

The transient solution of the system represented by equation (B-14a) 
can be written as follows: 

u' = 'e ^n*" sin ( /\-Cp- iv^t - cp) (B-70) 

From its definition in equation (B-28) , the logarithmic decrement 6 
can be determined as follows: 

6 - In e - = In e^“n T ■ T (B-71) 

e-O n <VT) 

where T is the period of damped oscillation, given by 

T = (2TTC/w n )(l-cV 1/2 (B-72) 


Thus , we now have 


6 = 2ttC (1-C 2 )' 1/2 


(B-73) 


For small damping (£^<<0. equation (B-73) can be replaced by the 
following approximate expression: 


I 6 % 2tt£ 


(B-74) 
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Combining equations (B-65) and (B-74), one obtains the following: 



(B-75) 


In view of equation (B-63) and since the specific energies are 
given by; 


U dv 3 V V ' U v iU/V> 


(B-76) 


we have 


C = U , / 4rrU 
dv s 


(B-77) 


However, the specific strain energy is 


U =ct 2 /2E 
v 


(B-78) 


A power- law-damping material is defined by equation (B-22). 
Combining this relation with equation (B-78), one obtains: 


C = (E C d /8T T )o m ' 2 


(B-79) 


Thus, it is clear that Q is independent of stress amplitude only for a 
KL (Kimball -Love 11) material (m=2). Then 


C = E C d /8TT 


or, in view of equation (B-65), 


(B-80) 


g = E C d /4rr 


(B-81) 


Furthermore, since the logarithmic decrement of a KL material with 
small damping is approximately proportional to the damping ratio as shown 
by equation (B-74), the following approximate relation holds: 


6 « E C d /4 


(B-82) 
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APPENDIX C 

DETAILS OF FORMULAS USED IN SECTION II 

Cl. Airy Stress Function and the Associated 
Stress, Strain, and Displacement Components 

Some useful formulas pertaining to analyses carried out in Section 
2.3 and 2.5 are summarized as follows. See the main text and the list 
of symbols for the definitions of notations used. It is to be noted that 
the same formula will apply to fiber as well as matrix regions with ap- 
propriate interpretations of the material property values, E and v, etc., 
and choices of terms in the series. For example, in the fiber region, 
terms containing negative powers of p must be omitted to prevent a 
singularity in stress. 

Polar stress components . - Stress components o ' , cJq, and are 
given by: 

-2 ' -5 

0 = a p + 2 b + 2 b. p cos 9 ■ 2 a.p cos 8 

r o o 1 1 

00 

- y [a n(n-l)p n 2 + b (n-2) (n+1) p n 

c_i n n 

n=2 ,3, . . . 

+ a^ n(n+l) p n 2 + b n (n+2)(n-l)p n ] cos n 9 (C-l) 
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-2 ' -3 

Og — a^p + 2 b Q + 6 b^p cos 9 + 2 a^p cos 9 

eo 

+ [a n n(n-l) p n ^ + b n (n+2) (n+1) p n +a^ n(n+l)p n ^ 

n=2 ,3, . . . 

+ b^ (n-2)(n-l) p n ] cos n 9 (C-2) 


' -3 

T r g = 2 bj p sin 9 - 2 a^p sin 9 

as 

+ ^ [a n n(n-l)p n ^ + b R n(n+l)p° 

n-2 ,3 , . . 

-a^ n(n+l) p n ^ - b n n(n-l)p n ] sin n 9 (C-3) 


Rectangular stress components . - Stress components o^.o^, and t 


are given by: 


-2 -3 

0 = a p cos 2 9 + 2 b + 2b,p cos 9 - 2 a ' p cos 3 9 

x o H o 1 K 1 K 

00 

-7 | a n(n-l)p n ^ cos (n-2) 9 + b (n+1) p n [n cos (n-2) 9 

t_i t n n 

n= 2 ,3 , . . . 

- 2 cos n 9] + a^ n(n+l)p n ^ cos (n+2) 9 

+ b'(n-l)p n [n cos (n+2) 9+2 cos n 9] ) (C-4) 

n j 

-2 -3 

a = -ap cos 2 0 + 2 b + 6 b.n cos 9 + 2 a! p cos 3 9 
y o p o 1 H 

00 

+ T (an (n-1) p° ^ cos (n-2) 9 +b (n+1) p n [n cos (n-2) 9 
Li L n n 

n= 2 , 3 , . . . 
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cos (n+2) Q 


+ 2 cos n 0] + a^ n(n+l)p 

> 

+ b^(n-l)p n [n cos (n+2) 0-2 cos n 0] | (C-5) 

T xy = a o P ^ sin 2 0 - 2 b^p sin 0 + 2 a^p ^ sin 0 (1-4 cos 2 0 ) 

00 

+ ^ [a R n(n-l)p n 2 + n(n+l)p n - a^ n(n+l) p n 2 
n=*2 ,3 , . . . 

+ b', n(n-l)p n ] sin (n+2) 0 (C-6) 

n 


Rectangular displacement components . - Displacement components u 


and v are given by: 


E (1+v) \u/R)=-a o p * cos 0 + 2 b Q (l-2v)p cos 0 

2 -2 
- b^p [3+4 (1+v)] cos 20+ a|p cos 2 0 

oo 

+ ^ {" a n npn ^ COS ® + ^n pn+ ^ [ -n COS ® 

n=2 , 3 , . . . 


+ 2 (l-2v) cos (n+1) 0-2 sin n 0 sin 0] 

+ a^ np n ^ cos (n+;l) 0 

+ b^p [n cos (n+1) 0+2 (l-2v) cos (n-1) 0 
+ 2 sin n 0 sin 0] ^ (C-7) 
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E (1+v) (v/R)- ~a o P ^ sin 9 + 2 b o (l-2v)p sin 9+4 b^p^(l-v) sin 2 9 

so 

+ ajp ^ sin 2 9+ ^ np n ^ sin (n-1) 9 

n=2 , 3 , . . . 

+ b^p n+ ^[ n sin (n-1) 9 + 2 (l-2v) sin (n+1) 9+2 sinn 9 cos 9] 

+ S n °P h ^ s * n ® 

+ b^p n+ ^ [n sin (n+1) 9 -2 (l-2v) sin (n-1) 9 -2 sin n 9 cos 9 ]| 

(C-8) 


C2 . Details of Formulas Used in Section 2.6 

Dirichlet torsion function Y*(o,9) (i-f.m) . - Torsion functions 
Y^ and y" 1 that satisfy Laplace equation in the fiber and matrix regions, 
respectively are assumed to be: 

00 

Y^(p,9) = a Q + ^ a k p k COS k ® 

k=l ,2 , . . . 

00 

Y m (p,9) - b o + ][ (b k p k + b_ k p" k ) cos k 9 

k=l 2 

Relationships among fiber-region and matrix-region coefficients a k 

and h k> - For the nth element, the distance of the fiber from the origin is 
d = 2(n-l)jj. The boundary conditions at the fiber-matrix interface are: 


(C-9) 


(C- 10) 
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XT f = V" + ^ (x-i) (% 2 W) 


dY f /dn = dV m /dn 


on C x } (G- 11) 


On 


the interface C^,p=l, hence, equations (C-li) becomes 


X(a o + Z cos k 0) = b^+^T b cos k ® 

k=l ,2 , . . . k=l ,2 , . . . 


+ (k) (X-l) (1+d + 2 d cos 9) 


I 

k=l ,2 , . . . 


k a. cos k 6 
k 


= ^ k (b^ - b cos k 0 
k=l ,2 , . . . 


(C-12) 

(C-13) 


Equating the coefficients of cos k 0 in equations (C-12 and 13) , we 
have : 


a = [b + (X-l)(l+d )/2]/X 
o o 


a 1 = 2 b 1 /(X+l) + X L d 


a k - 2 b k /(X+l) (k=2 ,3 , . . .) 


b_l= X 1 (b 1 -d) 


b -k = X l b k 


(k-2,3,...) 


(C-14) 


where 


X x a (X-1)/(X+1) 


(C-15) 
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Inter-element boundary conditions between the n-th and (n-t-l)-th 
elements . - On the inter-face boundary C 2 where | = (2 n-1) u , the 
local coordinate systems [p^ , 0^ .1)^ ] (i=l,2,...n) are inter- 

related as 


p (n)" p (n+l) * 9 (n) n ' 9 (n+l) ’ ^(n)" 11 ^(n) ?- 2 (n-l>U 

(C- 16) 


Thus the boundary conditions = Y™ n+ ^and = ^(n+1)^^ 

which warrant the displacement continuity and stress equilibrium on 
become, respectively: 


Cb, (n) -(-i)V n+1) Kol 


k=l , 2 .... 


= 2 X 1 (d-hi) 


p (i) 


cos 9 


(n) 


(C-17) 


I k[b k n) + ( - 1)kb k n+1)][p (n) COS (k ‘ l) 9 (n)- X l p (n)' 1 =os (k+1) 0 (n) ] 
k=l 2 

A. y ) • • • 


= 2\^u p cos 2 9^ 


(C-18) 


Shear stress formulas . - The shear stresses t *,t * (i=f,m) in 

x z yz 

the fiber and matrix regions of the n-th element are calculated from 


equations (79 and 80) (Section 2.6) as follows: 
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t f /(<*XG r) = OY f /BTD - Tl 
xz m 

00 

= - £ k a k o k_1 Sin (k-1) 6 -J\ (C-19) 

k=l ,2 .... 

T yf /(c ^ G m r) = " ^t f /S§) + l 

00 

= ' Z k a k pk_1 C0S (k_1) 9 + 5 (C-20) 

k=l ,2 , . . . 

T®/(oGr) = (av m /aT\) - T\ 
xz m 

00 

= “ Z k b k ^ pk_1 sin (k-1) 0 + x 1 p" k " lsin ( k+1 ) 0 1 

k=l ,2 , . . . 

+ (Xjd) p' 2 sin 2 0 - Tl (C-21) 

t “/(aGr) - - (&¥ m /35) + § 
yz m 

00 

= - Z k b k Cp k cos ^ - i) 9 ~ x^p k * cos ( k+ i) 0 3 

k=l ,2 , . . . 

+ (Xjd) p' 2 cos 2 0 + i (C-22) 


Torsional rigidity calculations for the nth element . - The torsional 
rigidity D ( n ) for the n-th element is calculated from equation (83 and 84) 
as 
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(C-23) 


D („) /(G m t > ' X < 2 VV + (2 W 


where 


R i ■ JT 't "5 OH - " ■„ 


R 2 = jj (| 2 +'h 2 ) d| dll = TT C (%) +d 2 ] 


R . JJ t" d 5 dH 


= b (4 6 - rr) + b 0 [(4/3) u 4 6 ( 1 - 6 2 ) + (n - arctan 6) ] 

o l X 


+ 4 £ b kM k+2 /(k+2) (Sj+ 6 k+2 S 2 ) 


k=4 , 6 , . . . 


+ 4 X 1 I b k u’ k+2 /(-k+2) (S 3 +6 _k+2 S A ) 


k=4 , 6 , . . 


R 4 = (4/3) u 4 6 (1+6 2 ) + 4u 2 6 d 2 - n [(%) + d 2 ] 


w. 


-1 

tan 6 


k+2 

sec 9 cos k 9 d 0 


C(k+2)/2] 


■ l 


(-w- 1 ** n ' 1 'ici)'* 2 -” 


V k+2 _ , „ 

S 2 = ) esc 9 cos k 0 d9 
tan"^6 
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C3. Details of Formulas Used in Section 2.7 

Saint-Venant flexural function x^(p» 0) (i=f,m).- Flexural functions 
and X m that satisfy the Laplace equation in the respective fiber and 
matrix regions are assumed in series form as follows: 


f V k 

X (P»0) = £ a k P cos k 9 

k=l , 3 

00 

X m (p,9) = Y, ^ b k pk + b -k p " k ^ C0S k 6 
k=l ,3 



Note that in equations (C-26) , coefficients with even-numbered 
subscripts are zero due to the anti-symmetry condition with respect to 
the T]-axis. 
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Relationships among fiber-region and matrix-region coefficients . 


a^ and b^. - In view of the displacement-continuity and the stress- 
equilibrium requirements at the fiber -matrix interface, where p = 1, 
a^ are related to b^, on substitution of eq. (C-26) into eqs. (109 and 
1 10) , as follows : 


a k = b k + b_ k (k = 1,3,...) 


Xa k = b k - b_ k (k = 5,7,...) 


Xa L = b x - b_ L - (X - 1) (3+2 v)/4 


Xa,=b - b + (X - l)/4 
J z -z 


> (C-27) 


Equation (C-27) may be solved for a k and b k in terms of b k to yield: 


aj = 2 b^/ (X+l) - X 1 (3+2v)/4 

a 3 = 2 b 3 /(X+l) + X x /4 

a k - 2 b k /(X+l) (k=3 , 5, . . .) 

b_ 3 = -Xj [b^ + (3+2v)/4] 

b_3 = -X L [b 3 + (1/4)] 

b_ k = - Xj^b k (k=3 , 5 , . . .) 

where X^ is as defined in eq. (C-15). 

Longitudinal thickness-shear stiffness S^. - In view of notations 
in eq. (106), the longitudinal thickness-shear stiffness expression, eq . 
(120), may be written as 


123 


S 55 /(E m r;> “ 4 6 M> 2 [ a'-D (tt/4) + 4 6 U 4 /3]/P' 


(C-29) 
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APPENDIX D 


EXPERIMENTAL OBSERVATION ON THE 
DAMPING BEHAVIOR OF A BORON FIBER 


In the damping analysis of a monofilament composite, it was found that 
the loss tangent associated with the in-plane Young's modulus is re- 

lated to the volume fractions, stiffness ratio, and the loss tangents of 
the constituent materials by the equation, [eq. (156)] 


‘11 


= (A 1 V g + V g ) /(A 1 v + V ) 
r E. m b t m 

f m 


In the case of boron-epoxy composite, the stiffness ratio is usually of 
the order of 120; whereas the loss tangent of boron is expected to be about 
one-tenth that of epoxy. The complete omission of the contribution of 
the boron fibers to the damping of the composite will then lead to a com- 
posite loss tangent which is unreasonably low. 

Unfortunately, so far as known to the authors, no experimental data 
dealing with the damping behavior of boron material alone are available 
in the literature. In view of this, a crude exploratory experiment was 
carried out on an Avco 4. 5-mil-diameter boron fiber in order to assess 
roughly the order of magnitude of damping in the boron material. A boron- 
fiber cantilever with a concentrated mass attached at its tip was deflected 
a certain prescribed distance and then released to oscillate in the vertical 


plane. The profile of the oscillation was recorded by a 16-mm motion-picture 
camera (set at 64 frames per second) until the motion of the fiber decays 
and returns to its initial equilibrium position. The decaying sinusoidal 
motion of the concentrated mass at the fiber tip is then reconstructed from 
the frame-by-frame observation of the film. The experiment was repeated 
several times with different initial deflections. 

The logarithmic decrement was obtained from the displacement versus 
time plots of the fiber tip (figure D-l) and then related to the loss 
tangent using eq. (B-75) from Appendix B. Averages of several runs were 
summarized as in figure D-2, where logarithmic decrements 6 based on the 
number of cycles elapsed were calculated for three cases with initial 
deflections equal to 1.5, 1.0, and 0.5 inch, respectively. These curves 
show clearly the amplitude dependency of the logarithmic decrement which 
is attributed mainly to the air damping. In view of the air damping which 
predominated, the loss tangent at small deflections, as estimated from 
this experiment, is of the order of 0.02 to 0.05, which is roughly ten 
times that of the estimated loss tangent for boron in vacuum (see Section IV) . 
It is concluded that, in order to assess the material damping of a boron 
fiber, the experiment must be carried out in vacuum to eliminate the 
effects of air damping which is both nonlinear and amplitude-dependent. 
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APPENDIX E 


COMPUTER PROGRAM DOCUMENTATION 


The computer program for computing the stiffnesses and the 
associated loss tangents consists of one lead-in program and ten 
subroutine programs. 

The lead-in program is concerned with the input of the material 
data, calling of each subprogram, and the output of the computed 
results . 

The input data consist of various pertinent material data 
such as Young's modulus, shear modulus, Poisson's ratio, and their 
respective loss tangents for specific frequencies as listed in Tables 

V-IX. 

Each subprogram is concerned with the calculation of a 
specific stiffness and associated loss tangent based on the input 
data. For example, subprogram C0MPE11 calculates the major Young's 
modulus and the associated loss tangent gg^. 

Each subprogram is designed to accomplish three functions: 

1. Generation of the parameters necessary for solution, 

2. Evaluation of the stiffness, and 

3. Evaluation of the loss tangent. 

The program was written in FORTRAN IV language as prescribed 


in IBM System Reference Library Form C-28-6274-3. 

A complete listing of the computer program is presented at 
the end of this report. 
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TABLE I. COMPARISON OF LONGITUDINAL IN-PLANE POISSON'S RATIO 


E,/E 

Longitudinal In-Plane Poisson's Ratio 

V f 

r tti 

Eq. (30) 

Eq. (38) 

120. 

0.331 

0.322 

0.3 


0.317 

0.311 

0.4 


0.316 

0.300 

0.5 


0.292 

0.288 

0.6 

5.3 

0.286 

0.282 

0.4 


0.273 

0.270 

0.5 


0.260 

0.258 

0.6 

24.1 

0.299 

0.294 

0.4 


0.284 

0.282 

0.5 


0.269 

0.268 

0.6 


5 


















TABLE II. LONGITUDINAL FLEXURAL STIFFNESS EFFICIENCY OF 


MONOFILAMENT COMPOSITES. 


Composite 

E c /E 
f m 

■ 6 

U 

V f 


E (b) /E 
E 11 /E 1I 


Eq. (48) | Eq. (142) 


Steel-Epoxy 

74 

1.00 

1.11 

0.636 

0.616 



1.00 

1.06 

0.708 

0.683 



0.95 

1.11 

0.673 

0.683 



0.90 

1.11 

0.707 

0.755 



0.95 

1.06 

0.746 

0.755 



0.93 

1.16 

0.743 

0.730 



1.00 

1.00 

■ 0.785 

0.755 

S glass-Epoxy 

24 

0.85 

1.38 

0.482 

0.578 



0.90 

1.23 

0.572 

0.636 



0.95 

1.11 

0.673 

0.697 

Boron-Epoxy 

120 

0.85 

1.38 

0.482 

0.549 



0.90 

1.23 

0.572 

0.614 



0.95 

1.11 

0.673 

0.681 

Boron-Al . 

6 

0.85 

1.38 

0.482 

0.678 


0.608 

0.677 

0.677 

0.750 

0.750 

0.725 

0.750 


0.542 

0.608 

0.677 



6 




























TABLE III. COMPARISONS OF POISSON FLEXURAL STIFFNESS 


FOR COMPOSITES HAVING SQUARE TYPICAL ELEMENTS. 


Composite 

Fiber 
Volume 
Fract ion 

^ D 12^eq. (50) / ^ D 12^eq. (4) 

Boron-Epoxy 

0.4 

0.725 


0.5 

0.593 


0.6 

0.567 

Glass-Epoxy 

0.4 

0.780 


0.5 

0.657 


0.6 

0.638 

Bor on -A 1 . 

0.4 

0.930 


0.5 

0.835 


0.6 


0.833 












TABLE IV. COMPARISONS OF IN-PLANE LONGITUDINAL SHEAR MODULUS AND 


EQUIVALENT SHEAR MODULUS AS CALCULATED FROM THE TORSION ANALYSIS. 


Composite 


G./G 
f m 


Fiber 

volume 

Ratio G,,/G 
66 m 

fraction 

In-Plane 

Equivalent in Torsion 

0.4 

2.30 

25.30 

0.5 

3.23 

38.98 

0.6 

4.67 

55.66 

0.4 

2.16 

5.04 

0.5 

2.84 

7.34 

0.6 

3.80 

. 

10.13 

0.4 

1.84 

2.01 

0.5 

2.22 

2.58 

0.6 

2.72 

3.29 


Boron-Epoxy 135. 


E Glass-Epoxy 23.4 
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TABLE V, - ELASTIC AND DAMPING PROPERTIES OF BORON 
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TABLE VI. - ELASTIC AND DAMPING PROPERTIES OF EPOXY A 
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TABLE VII. - ELASTIC AND DAMPING PROPERTIES OF EPOXY B 
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TABLE VIII. - ELASTIC AND DAMPING PROPERTIES OF ALU 2224-T3 
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TABLE IX. - ELASTIC AND DAMPING PROPERTIES OF E-GLASS 
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TABLE X. - IN-PLANE LONGITUDINAL YOUNG'S MODULUS AND ASSOCIATED LOSS TANGENT 
FOR A BORON-ALUMINUM 2024-T3 MONOFILAMENT COMPOSITE 
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TABLE X. - IN-PLANE LONGITUDINAL YOUNG’S MODULUS AND ASSOCIATED LOSS TANGENT 
FOR A BORON-ALUMINUM 2024-T3 MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE X. - IN-PLANE LONGITUDINAL YOUNG'S MODULUS AMD ASSOCIATED LOSS TANGENT 
FOR A BORON-ALUMINUM 2Q24-T3 MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XII. - IN-PLANE TRANSVERSE YOUNG'S MODULUS AND ASSOCIATED LOSS TANGENT FOR 
A BORON -ALUMINUM 2024-T3 MONOFILAMENT COMPOSITE 
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TABLE XII. - IN-PLANE TRANSVERSE YOUNG 'S MODULUS AND ASSOCIATED LOSS TANGENT FOR 
A BORON -ALUMINUM 2024-T3 MONOFILAMENT COMPOSITE (CONTINUED) 
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TABL- XII - IN-PLANE TRANSVERSE YOUNG'S MODULUS AND ASSOCIATED LOSS TANGENT FOR 
A BORON-ALUMINUM 2024-T3 MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XIII. - IN-PLANE LONGITUDINAL SHEAR MODULUS AND ASSOCIATED LOSS TANGENT FOR 
A BORON-ALUMTNIJM 202^-T3 MONOFILAMENT COMPOSITE 
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TABLE XIV. - LONGITUDINAL FLEXURAL STIFFNESS AND ASSOCIATED LOSS TANGENT FOR 
A BORON -ALUMINUM 2024-T3 MONOFILAMENT COMPOSITE 


X oooooooooooo 

** UIUJUJUJUIUJUJUIUJUJUJUJ 

o oo«fff«»N\otn«- 

-4 '0'C®tnininwuiininintn 

o «••••••••••• 

oooooooooooo 


o 




o 




4 




• 




o 






X 


u. 


UJ 


> 


V 




u. 




tu 



• 




• 



CJ 

* 



o 

* 



o 

* 



• 





►* 




t 




4 




CM 



< 

O 



►- 

CM 


z 

-1 


wm 

Ul 

111 

X 

. •« 

u 

a 

-J 

a 



z 

. o 

a 




Ul 


X 


a 


3 




_) 



— • 

< 



o 

1 



4 

z 



• 

o 



—4 

a 




a 




ai 




« 




* 


H 

X 

» 


tn 

a 

* 

in 

a 


* 

w 




44 

>0 



o 

Ul 




n 


44 

o 

o 

•4 

o 

44 

o 

** 

o 

44 

O 

•4 

o 

O 

44 

o 

o 

o 

O 

Ul 

Ul 

Ul 

Ui 

UJ 

Ul 

Ul 

Ul 

Ui 

Ul 

Ul 

Ul 

CM 

a* 

o 

® 

CM 

Ul 


*4 

<1 

in 


N 

IT 

CM 

44 

a 

CM 

CO 

Ul 

44 

N 

vO 

in 

m 

CM 


o 

® 

N 

•o 


<o 

in 

Ul 

in 

Ul 

IT 

i n 

Ul 

4 

4 

4 

4 

♦ 

4 

4 


4 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 


o 

O 

o 

o 

o 

o 

o 

o 

*4 

m4 

44 

.4 

o 

o 

o 

o 

o 

o 

o 

a 

o 

o 

a 

u 

Ul 

Ul 

Ul 

Ui 

Ul 

Ul 

Ul 

Ul 

1 

UJ 

1 

UJ 

9 

UJ 

1 

Ui 

44 

1*1 

o 

o 


CM 

CM 

4 

N 

N 


® 

n 

Ul 

*4 

4 

r» 

O 

4 


CM 

CM 

n 

•4 

44 

CM 

4 

N 

in 

4 

CM 

CM 

® 

N 

>o 

N 


*4 

44 

44 




* 

CD 

CD 

® 

® 

O 

o 

o 

O 

o 

o 

O 

O 

O 

O 

o 

O 


CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

cm 

CM 

eg 

CM 

O 

O 

O 

O 

O 

O 

O 

O 

o 

O 

o 

O 

Ul 

UJ 

Ul 

tu 

Ul 

Ul 

Ul 

UJ 

UJ 

Ul 

UJ 

UJ 

<0 

4 

f" 

<c 

in 

O' 

< 

«o 

o 

O' 

N 

4 

to 

4 

CVI 

o 

O' 

® 

® 

K 

s 



>o 

N 

N 

*4 

r* 

44 

b- 

*> 

44 

<£> 

vO 

44 

44 


U) 

44 

MD 

• 

O 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 


< 

Q 




44 

•» 

CM 

CM 

<M 

cm 

"> 

n 

n 

ro 

** 


> 


O 

o 

O 

O 

O 

C 

o 

o 

o 

o 

o 

o 

U 














z 

N 

Ul 

UJ 

UJ 

UJ 

JJ 

Ui 

IU 

ui 

UJ 

UJ 

UJ 

Uj 

JJ 

►- 

O 

o 

o 

o 

o 

O 

o 

c 

c 

o 

o 

o 

J 

a 

o 

o 

o 

o 

o 

o 

3 

o 

c 

o 

o 

o 

<J 

UJ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

UJ 

i 

in 

N 

— « 

CM 

in 

K 

44 

CVJ 

in 


*4 

CM 

a 


• 

• 

• 

• 

• 

• 

• 

4 

• 

• 

• 

• 

u. 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 


166 



TABLE Xiv. - LONGITUDINAL FLEXURAL STIFFNESS AMD ASSOCIATED LOSS TANGENT FOR 
A BORON -ALUMINUM 2024-T3 MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XIV. - LONGITUDINAL FLEXURAL STIFFNESS AND ASSOCIATED LOSS TANGENT FOR 
A BORON -ALUMINUM 20?A-T3 MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XV - POISSON FLEXURAL STIFFNESS AND ASSOCIATED LOSS TANGENT FOR 
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TABLE XVI. - TRANSVERSE FLEXURAL STIFFNESS AND ASSOCIATED LOSS TANGENT FOR 
A BORON -ALUMINUM 2024-T3 MONOFILAMENT COMPOSITE 
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TABLE XVI. - TRANSVERSE FLEXURAI. STIFFNESS AND ASSOCIATED LOSS TANGENT FCR 
A BORON -A LIT'TNUM 2024-T3 MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE! XVII. - TWISTING STIFFNESS AND ASSOCIATED LOSS TANGENT FOR A BORON - 
ALUMINUM 2024-13 MONOFILAMENT COMPOSITE 
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TABLE XVIT - TWISTING STIFFNESS AND ASSOCIATED LOSS TANGENT FOR A BORON- 
ALUMINUM 2024-T3 MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XVII. - TWISTING STIFFNESS AND ASSOCIATED LOSS TANGENT FOR A BORON - 
ALUMINUM 2024-T3 MONOFILAMENT COMPOSITE (CONTINUED) 
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table xvm. - longitudinal thickness^shear modulus and associated loss tangent 

FOR A BORON-ALl'MINUM 202H-T3 MONOFILAMENT COMPOSITE 
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TABLE XVIII. - LONGITUDINAL THICKNESS -SHEAR MODULUS AND ASSOCIATED LOSS TANGENT 
FOR A BOP.ON -ALUMINUM 2024-T3 MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XIX. - TRANSVERSE THICKNESS-SHEAR MODULUS AND ASSOCIATED LOSS TANGENT FOR 
A BORON-ALUMINUM 2024-T3 MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XX. - IN-PLANE LONGITUDINAL YOUNG'S MODULUS AND ASSOCIATED LOSS TANGEN 1 
FOR AN E GLASS-EPOXY MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XXT. - IN-PLANE MAJOR POISSON'S RATIO AND ASSOCIATED LOSS TANGENT FOR 
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TABLE XXIII. - IN-PLANE LONGITUDINAL SHEAR MODULUS AND ASSOCIATED' LOSS TANGENT FOR 
AN E GLASS -EPOXY MONOFILAMENT COMPOSITE 
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TABLE XXIII. - IN-PLANE LONGITUDINAL SHEAR MODULUS AMD ASSOCIATED LOSS TANGENT FOR 
AN E GLASS-EPOXY MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XXIII, - IN -PLANE LCNGITU OINA L SHEAH MODULUS AND ASSOCIATED LOSS TANGENT FOR 
AN E GLASS-EPOXY MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XXIV. - LONGITUDINAL FLEXURAL STIFFNESS AND ASSOCIATED LOSS TANGENT FOR 
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TABLE XXIV. - LONGITUDINAL FLEXURAL STIFFNESS AND ASSOCIATED LOSS TANGENT FOR 
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TABLE XXIV. - LONGITUDINAL FLEXURAL STIFFNESS AND ASSOCIATED LOSS TANGENT FOR 
AN E GLASS-EPOXY MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XXVI. - TRANSVERSE FLEXURAL STIFFNESS AND ASSOCIATED LOSS TANGENT FOR 
AN E GLASS-EPOXY MONOFILAMENT COMPOSITE 
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TABLE XXVT. - TRANSVERSE FLEXURAL STIFFNESS AND ASSOCIATED LOSS TANGENT FOR 
AN E GLASS-EPOXY MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XXVI. - TRANSVERSE FLEXURAL STIFFNESS AND ASSOCIATED LOSS TANGENT FOR 
AN E GLASS-EPOXY MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XXVII. - TWISTING STIFFNESS AND ASSOCIATED LOSS TANGENT FOR AN E GLASS- 
EPOXY MONOFILAMENT COMPOSITE 
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TABLE XXVII. - TWISTING STIFFNESS AND ASSOCIATED LOSS TANGENT FOR AN E GLASS- 
EPOXY MONOFILAMENT COMPOSITE (CONTINUED) 
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TABLE XXVII. - TWISTING STIFFNESS AND ASSOCIATED LOSS TANGENT FOR AN E GLASS - 
EPOXY MONOFILAMENT COMPOSITE (CONTINUED) 
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table xxix. - transverse thickness -shear modulus and associated loss tangent for 

AN E GLASS-EPOXY MONOFILAMENT COMPOSITE 
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(a) Pure bending of a typical element 





figure 5. Longitudinal flexural loading of a typical mono- 
filament composite element. 
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(a) Tip loading of a monofilament beam 



5 


(b) Cross-sectional view 

Figure 8. Tip loading of a monofilament composite 


element , 









IMAGINARY AXIS 



Figure 10. Concept of complex moduli end 
loss angle. 
















Figure 13. Transverse flexural stiffness 
efficiency. 


225 




Figure 14. Comparisons of Prandtl torsion 
function (0x10^) of a square 
matrix with a circular insert. 

Remarks: Numbers at the mesh 

points denote the values of 
0xl(A obtained by Ely and 
7 ienkiewics (ref. 75) and 
those from the point -matching 
method (in parentheses). 
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ELEMENT— +• FIRST ELEMENT 



three-element model 










Figure 18. Transverse thickness -shear modulus. 
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Figure 19. Dynamic Young's modulus and associated loss tangent for boron 



(Curve A) 



Figure 20. Dynamic Young's modulus and associated loss tangent for epoxy 
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(Curve A) 



sson's ratio and associated loss tangent 
monofilament composite. 
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-plane transverse Young's modulus and associated loss 
ngent far a boron-epoxy monofilament composites. 






Q 

CM 


99 ^ 
(%) ^ 


1N39NV1 SSOT 



a 


a 






(isd^oo"9 smnacw uvshs ivmmmoi 3N\nd-Ni 



a 
d 
nj 

r-i C 

CL O 


I 

c 


o 


M U3 


r-. 

<nj 

<U 

D 

00 

•r^ 

Cxh 


239 


-epoxy monofilament composite. 
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-epoxy monofilament composite. 
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epoxy monofilament composite. 







Figure 31. Twisting stiffness and associated loss tanget 
for a boron-epoxy monofilament composite. 
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Figure 32. Longitudinal thickness-shear stiffness and associated 
- - loss tangent for a boron-epoxy monofilament composite 
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boron-epoxy monofilament composite. 
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THREE-PARAMETER MODELS 
(STANDARD LINEAR SOLID) 


Figure B- 1 . Two-and three- parameter viscoelastic models. 
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